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SUMMARY 


The  primary  objective  of  this  project  is  to  develop  and  apply  improved  statistical  methodologies  for 
nuclear  monitoring.  During  the  past  two  years,  our  major  efforts  have  been 

[A]  Measure  the  teleseismic  P- wave  magnitudes  of  historical  Soviet  explosions  as  well  as  explosions 
from  other  foreign  test  sites  recorded  at  the  optimal  distance  range  from  20°  to  95°. 

[B]  Perform  various  statistical  analyses  of  the  raw  mb  and  obtain  the  optimal  network  mb  values.  Con¬ 
duct  the  maximum-likelihood  magnitude-yield  regressions  and  analyze  the  source-depth  scaling  relation¬ 
ship. 

[C]  Conduct  a  theoretical  study  to  investigate  relevant  issues.  Improve  and  document  the  statistical  as 
well  as  the  forward-modeling  tools  currently  in  use. 

This  semi-annual  report  summarizes  our  results  obtained  under  Task  [C]  during  the  period  Sep¬ 
tember  20, 1993,  through  November  5, 1993. 

An  accurate  prediction  of  the  regional  phases  due  to  explosion  sources  in  areas  of  high  prolifera¬ 
tion  concern  requires  a  decent  understanding  of  the  attenuatiorVscattering  mechanisms  along  the  propa¬ 
gation  paths.  Motivated  by  many  observations  of  correlation  between  the  free-surface  topography  and 
the  coda  generation  at  teleseismic  and  regional  distances,  seven  simple  statistics  have  been  presented 
and  analyzed  in  this  study  to  test  whether  they  can  fully  characterize  the  roughness  of  an  arbitrary 
rugged  free-surface  topography.  The  definitions  of  these  parameters  are  very  natural,  and  some  of 
parameters  might  be  suitable  to  quantify  the  roughness  of  realistic  profiles.  The  normalized  arc  length 
and  Durbin-Watson  number  (von  Neumann  ratio)  are  identified  as  the  best  and  the  poorest,  respectively, 
roughness  measures,  in  terms  of  correlating  with  the  propagation  behavior  of  Rg  (see  Figures  95-102 
for  the  regressions).  2-dimensional  finite-diflerence  simulations  are  then  conducted  to  investigate  short- 
period  Rayleigh  wave  scattering  by  these  27  topographic  profiles.  The  numerical  experiments  show  that 
rougher  topography  causes  stronger  attenuation  (as  expected),  that  reflection  is  an  efficient  process, 
and  that  the  bulk  of  the  energy  which  is  not  transmitted  across  the  topography  is  either  converted  to 
body  waves  or  reflected  back.  The  scattered  P  and  S  waves  are  radiated  as  coda  to  teleseismic  dis¬ 
tances  and  some  P-SV  energy  is  trapped  in  the  crustal  waveguide  as  Pg ,  Sg,  and  Lg .  The  body 
waves  appear  to  radiate  from  the  changes  in  slope  and  to  first  order  can  be  visualized  as  from  point 
diffractors  along  the  rugged  free-surface.  For  most  models  tested  in  this  study,  the  Rg  -to -S  conversion 
is  more  efficient  than  that  of  Rg  -to-P  (see  Figures  103-129).  Particle  motion  plots  of  the  scattered 
bodywaves  can  help  to  illustrate  this  point  too  (see  Figure  130-156).  The  displacement  synthetics  of  the 
scattered  bodywaves,  which  makes  the  synthetics  can  be  compared  easily.  This  approach  is  more  con¬ 
venient  than  another  one  which  presents  the  scattered  wavefield  in  strain  units.  The  spatial  O  values 
associated  with  rough  topographic  profiles,  inferred  from  the  logarithm  of  the  transmitted  power  spec¬ 
trum,  range  from  15  to  80,  which  are  in  excellent  agreement  with  the  Op  of  the  uppermost  crust  in  many 
places  of  the  world.  Regression  results  also  suggest  that  both  the  attenuation  coefficient,  y,  and  the 
scattering  loss  can  be  represented  as  linear  functions  of  the  roughness  parameters.  This  work  offers  a 
promising  approach  to  quantify  effects  of  scattering  by  topography  and  heterogeneity  upon  other 
regional  phases. 
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1.  INTRODUCTION 


The  problem  of  short-period  fundamental-mode  Rayleigh  waves  (or  Rg )  scattered  by  topography  is 
important  and  interesting  for  many  reasons.  A  significant  fraction  of  the  elastic  energy  generated  by  a 
shallow  explosion  is  contained  in  short-period  Rayleigh  waves.  These  short-period  Rayleigh  waves  are 
the  largest  phase  on  the  near-source  seismogram.  However,  very  often  they  are  attenuated  rather 
rapidly  with  spatial  (7s  at  1  Hz  in  the  range  of  20  to  200.  Such  attenuation  can  be  due  to  intrinsic  shear 
attenuation  in  the  near  surface,  or  it  can  be  accomplished  by  scattering  from  shallow  heterogeneity  and 
topography.  If  the  energy  is  scattered  then  it  must  be  observed  as  coda  in  the  near-source,  regional, 
and/or  teleseismic  seismogram. 

The  scattering  of  short-period  Rayleigh  waves  by  topography  has  been  proposed  as  a  mechanism 
for  the  generation  of  teleseismic  P  coda  by  Greenfield  (1971).  Gupta  and  Blandford  (1987)  (and  Gupta 
and  Wagner,  1992)  propose  that  such  P  coda  may  be  diagnostic  of  the  near  source  medium.  The  abil¬ 
ity  of  using  teleseismic  complexity  for  discriminating  explosions  from  earthquakes  can  be  accounted  for 
by  an  explanation  that  the  emergent  earthquake  teleseismic  coda  is  comprised  of  pP  and  of  S-to-P 
conversions  at  relatively  sharp  interfaces,  whereas  the  emergent  explosion  coda  is  comprised  of  Rg  -to- 
P  conversions.  Thus  the  coda  in  emergent  explosion  signals  would  have  more  of  a  “rolling"  character 
which  can  be  distinguished  from  the  “impulsive"  coda  of  emergent  earthquakes  (Blandford,  1993). 
Adushkin  et  al.  (1990)  find  Rg  from  Soviet's  Sernpalatirtsk  nuclear  test  site  a  very  stable  yield  estima¬ 
tor.  With  the  focus  of  nuclear  monitoring  research  moving  from  the  yield  estimation  to  seismic  event 
identification,  an  assessment  needs  to  be  made  of  the  maximum  range  at  which  Rg  could  be  used  to 
discriminate  small  events  in  various  geologic/tectonic  environments  (DARPA,  1992).  P  coda  observed 
at  arrays  is  often  partly  composed  of  Rayleigh  waves  generated  by  P  waves  steeply  incident  upon  shal¬ 
low  heterogeneity  or  topography.  Key  (1967  and  1968)  identified  such  Rayleigh  waves  within  the  P 
coda  at  the  EKA  array  (Eskdalemuir,  Scotland)  for  certain  azimuths  of  arrival.  Dainty  (1986)  shows  by 
use  of  frequency-wavenumber  analysis  that  P  coda  at  the  small  NORESS  array  (southern  Norway)  is 
often  composed  of  fundamental  short-period  Rayleigh  waves.  Both  the  coherency  of  such  coda  and  its 
persistence  are  closely  related  to  the  scattering  of  the  Rayleigh  waves  by  topography  and  shallow 
heterogeneity.  Aki  and  Chouet  (1975)  proposed  three  simple  models  for  the  generation  and  propagation 
of  short-period  coda  from  local  earthquakes.  They  considered  the  case  of  scattered  surface  waves  as 
well  as  scattered  body  waves.  Although  the  coda  of  local  earthquakes  appears  to  be  dominated  by  S 
waves,  it  is  clear  that  trapped  modes  and  surface  waves  are  also  present  (Phillips  and  Aki,  1986).  A 
portion  of  the  regional  Lg  phase  and  Its  coda  is  composed  of  fundamental  and  higher  mode  Rayleigh 
waves  trapped  in  the  crustal  waveguide.  Der  et  al.  (1984)  propose  that  a  significant  portion  of  the  Lg 
coda  is  composed  of  short-period  fundamental  Rayleigh  waves  that  concentrate  energy  near  the  surface 
and  are  multiply  scattered,  contributing  to  the  long  Lg  coda  observed  on  regional  seismograms. 
Baumgardt  (1990)  proposes  that  Lg  from  Eastern  Kazakh  incident  upon  the  Urals  produces  observable 
teleseismic  P  coda  at  NORSAR.  All  of  these  problems  have  in  common  a  need  to  estimate  the  effects 
of  scattering  on  short-period  surface  waves  due  to  shallow  heterogeneity  and  topography. 


Topographic  scattering  through  a  right-angle  step  was  addressed  by  deBremaecker  (1958)  and 
continues  as  a  topic  of  interest.  Numerous  workers  have  examined  the  incident  Rayleigh  or  Love  wave 
upon  welded  quarter  spaces  (see  for  example,  Mai  and  Knopoff,  1965;  McGarr  and  Aisop,  1967;  Drake, 
1972;  Munasinghe  and  Famell,  1973;  Tuan  and  Li,  1974;  Martel  et  al.,  1977;  Sabina  and  Willis.  1977; 
Chen  and  Aisop,  1979;  Fuyoki  and  Matsumoto,  1980;  Drake  and  Bolt,  1980;  Momoi,  1985;  Fuyuki  and 
Nakano,  1984;  Wojdk  and  Mak,  1985;  and  Nakano  et  al. ,  1988).  Ultrasonic  modeling  was  used  by  Bul- 
lit  and  Toksoz  (1985)  and  Toksoz  et  al.  (1986)  to  examine  the  obliquely  incident  Rayleigh  wave 
transmission  and  reflection  by  sloping  topography. 

Most  previous  research  concerned  with  the  scattering  of  Rayleigh  waves  has  dealt  with  the  sim¬ 
plest  geometric  configurations  (quarter-space,  three-quarter  space,  a  trench  or  a  vertical  step)  where  the 
boundaries  are  always  parallel  or  perpendicular  to  the  grid  axis.  The  present  study  concentrates  on 
more  general  topographic  profiles,  as  a  follow-up  of  our  previous  study  summarized  in  McLaughlin  and 
Jih  (1986  and  1987).  Linear  Finite-Difference  (LFD)  modeling  is  chosen  simply  out  of  convenience.  The 
method  has  the  advantages  that  the  solution  contains  all  conversions  (P-SV,  SV-P)  and  all  orders  of 
multiple  scattering.  We  use  the  2nd-ord er  formulation  of  Kelley  et  al.  (1976)  with  absorbing  boundary 
conditions  described  by  Clayton  and  Engquist  (1977)  and  Ist-order  free-surface  boundary  conditions  for 
topographic  surfaces  described  by  Jih  et  al.  (1988)  and  Jih  (1993a).  The  method  permits  examination  of 
fairly  general  2-D  models.  Surface  waves  are  introduced  as  initial  conditions  within  the  half-space  grid. 
There  is  now  a  considerable  body  of  work  covering  many  aspects  of  scattering  using  2-D  LFD  methods. 

McLaughlin  and  Jih  (1987)  examined  the  interaction  of  incident  Rayleigh  waves  upon  a  2-D  shal¬ 
low  heterogeneous  layer  of  various  thicknesses,  scale  lengths,  and  rms  velocity  fluctuations  embedded 
in  a  homogeneous  half  space  using  2-dimensional  LFD.  In  this  paper,  we  explore  the  stability  of  planar 
Rayleigh  waves  incident  upon  2-D  topography  using  the  same  method  of  analysis.  The  goal  is  to  deter¬ 
mine  some  basic  quantities  such  as  scattering  strength,  conversion  efficiency,  and  backscattering 
strength  of  the  topographic  relief.  These  quantities  are  needed  in  order  to  compose  models  for  coda 
observed  near  the  surface  of  the  Earth.  Although  the  present  study  tackles  the  same  problem  as  was 
treated  in  our  previous  work  (McLaughlin  and  Jih;  1986)  (hereafter  [1]),  there  are  several  fundamental 
improvements.  No  clear  (i.e. ,  quantitative)  relationship  between  the  roughness  measures  and  the  pro¬ 
pagation  characteristics  was  found  in  [1],  because  the  rough  topographic  profiles  used  in  [1]  did  not 
cover  a  wide  range  of  profile  types.  The  scattered  body  waves  were  described  in  non-directional  strain 
in  [1],  which  made  the  synthetics  more  difficult  for  direct  comparison.  In  the  present  work,  more  rough¬ 
ness  measures  are  tested  with  a  suite  of  realistic  profiles.  Furthermore,  the  scattered  body  waves  are 
presented  in  the  same  displacement  unit,  with  direction  of  each  ray  clearly  identified. 

We  first  examine  the  problem  of  transmission  and  the  partitioning  of  energy  by  the  scattering  pro¬ 
cess.  It  is  demonstrated  as  expected  that  rougher  topography  produces  larger  attenuation.  Reflection 
(backscattering)  is  observed  to  be  an  efficient  process  for  all  cases  examined.  In  this  way  the  rough 
topography  is  characterized  by  a  spatial  O  and  several  other  roughness  statistics.  We  use  the  con¬ 
verted  body  wave  field  at  depth  in  the  half-space  as  a  measure  of  the  simulated  far-field  P  and  S  coda. 
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Attenuation  of  the  short-period  Rayleigh  wave  is  found  to  be  strong  by  moderately  rough  topography, 
and  Rayteigh-to-body  wave  conversion  is  seen  to  dominate  the  attenuation  process.  Examination  of  the 
converted  body  wave  field  in  the  half-space  leads  us  to  conclude  that  a  significant  fraction  of  the  con¬ 
verted  SV  energy  will  be  trapped  in  the  crust  as  Lg . 


2.  ROUGHNESS  STATISTICS 


2.1  Generation  of  Random  Topography 

In  this  study,  it  is  assumed  that  the  rough  topography  can  be  approximated  by  level  elevations 
separated  by  sections  of  45°  slope.  Topographic  profiles  are  generated  with  a  random  number  generator 
which  produces  a  Markov  sequence  of  values  0,  1,  and  -1,  according  to  an  arbitrarily  assigned  probabil¬ 
ity  of  the  “state"  1,  hereafter  denoted  as  Pr|l].  We  assume  0  <  Pr[-1 1  =  Pr(1]  <  0.5  and  Pr(0)  -  1  -  Pr(1] 
-  Prf-l].  The  sequence  of  0, 1  and  -1  is  regarded  as  the  successive  elevation  change  in  the  topography. 
That  is,  the  integral  of  this  Markov  sequence  represents  a  profile  suitable  for  forward  modeling  experi¬ 
ments.  Larger  Pr[1]  would  cause  more  frequent  elevation  changes  in  the  resulting  topographic  profile, 
and  hence  it  is  an  indication  of  rougher  topography,  in  this  study,  both  smooth  and  extremely  rugged 
topography  are  examined,  with  Pr(1]  ranging  from  5%  to  35%.  Initially,  the  values  1  and  -1  are 
assumed  to  have  the  same  probability.  Nevertheless,  the  selection  of  the  current  “state”  of  the  Markov 
chain  is  governed  not  only  by  the  preset  probability  but  also  by  the  previous  state.  For  instance,  given  a 
previous  value  of  1 ,  the  current  value  can  only  be  0  or  1 .  This  Markov  condition  is  imposed  to  make 
sure  that  the  topography,  as  a  piecewise-linear  polygon,  has  at  least  two  subsegments  of  the  same 
slope  in  each  of  its  linear  segment,  and  to  prevent  the  occurrence  of  sudden  changes  such  as  from  1  to 
-1  and  vice  versa. 

27  rough  topographic  profiles  thus  generated  were  inserted  into  the  grid  with  flat  ground  con¬ 
catenated  on  both  sides.  Figure  1  through  9  show  the  rugged  profile,  H(x).  of  the  27  models  on  the  left, 
log-amplitude  wavenumber  spectra  of  the  topography  in  the  center,  and  several  corresponding  rough¬ 
ness  statistics  (which  are  to  be  defined  in  the  following  section)  listed  on  the  right  of  topographic  profiles 
numbered  11  through  73.  Only  the  nigged  portion  of  each  profile  is  shown.  In  all  these  models,  the 
rugged  portion  has  a  lateral  span,  X,  approximately  30  km. 

2.2  Definition  of  Various  Roughness  Measures 

In  searching  for  quantitative  measures  to  characterize  the  roughness  of  a  topographic  profile,  one 
would  expect  to  seek  parameters  which  are  easy  to  compute  once  the  profile  is  digitized.  Note  that  each 
discretized  profile  is  just  a  time  series  of  equal  spatial  spacing  in  the  lateral  direction,  and  therefore 
many  techniques  used  in  time  series  analysis  might  be  suitable  for  our  need.  Several  statistics  which 
appear  to  be  natural  candidates  of  the  roughness  measure  are  defined  and  tested  in  this  study. 

The  " normalized  total  variation”,  £2,,  is  defined  as  £2,  =  ^lHrHj_1l/X  where  Hj  is  the  elevation  (in 
km)  at  the  ith  column  of  the  nontrivial  portion  of  the  topography.  This  statistic  reflects  the  idea  that  in  a 
smooth  time  series  the  differences  between  successive  points  are  all  small,  and  hence  the  numerator 
will  be  small.  The  denominator  is  simply  a  scaling  factor.  £2,  is  nothing  but  the  "normalized  L1  norm”  of 
the  elevation  change.  This  parameter  turns  out  to  be  equivalent  to  the  total  "probability  of  elevation 
change"  of  the  topographic  profile.  Another  roughness  measure,  which  is  very  similar  to  £2,.  is  the  "nor¬ 
malized  L2  norm”  of  the  elevation  change,  £22,  defined  as  VS(H~Hi-i)2/X.  When  each  elevation  change 
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is  limited  to  have  at  most  one  unit  step  (as  in  all  the  examples  used  in  this  study).  ft2  is  equivalent  to  ft, 
(up  to  a  constant).  This  is  not  true  in  the  more  general  case,  and  hence  both  quantities  ought  to  be 
tested.  The  " Root-mean-squared  elevation",  ft^  is  defined  as  VZ(  H~H  )2 1  N  H  the  average 
elevation  of  the  topography,  and  N  the  number  of  subintervals.  Q4  e  £(H-Hj_,)2/£(Hr-H)2,  is  also  called 
the  “Durbin-Watson  statistic"  or  the  “ von  Neumann  ratio"  (cf.  Bloomfield,  1975).  The  peak  wavenumber 
spectrum,  K(0),  is  denoted  as  fts.  The  " number  of  slope  changes  per  unit  length"  (i.e. ,  the  average 
number  of  corners  per  km)  is  denoted  as  ft*  Note  that  ft«  is  directly  proportional  to  the  "probability  of 
slope  change",  whereas  ft,  reflects  the  "probability  of  elevation  change".  In  calculus  terms,  ft,  and  fte 
are  closely  related  to  the  first  and  second  spatial  derivatives,  respectively,  of  the  topographic  profile,  H. 
ft5,  the  peak  wavenumber  spectrum  K(0),  is  a  measure  of  the  mean  elevation,  H,  whereas  ft3  is  a 
measure  of  the  deviation  from  the  mean.  The  arc  length  (normalized  by  the  total  span),  ft7,  is  computed 
as  £V[ AX2  +  (HrHj_,)2]/X  where  AX  is  the  spacing  (in  km)  in  the  horizontal  direction.  For  a  flat  profile, 
ft7  equals  to  1.  There  exist  simple  relationships  among  some  of  the  roughness  measures.  For  instance, 

ft2  =  VA2  •  ft,/X  ;  and  ft4  -  N  •  (AX  •  ft2  /  ft3)2. 

Table  1  lists  various  roughness  measures  of  27  profiles  numbered  from  11  to  73  with  increasing 
Pr(1].  In  each  case,  the  random  number  generator  is  initialized  with  the  seed  number  being  the  trun¬ 
cated  system  clock  time  at  the  time  the  profile  was  generated.  These  seed  numbers  are  also  included 
in  Table  1,  and  hence  all  the  profiles  can  be  duplicated  on  SUN4  computers.  The  preset  probability  of 
state  1  varies  from  5%  (models  11  through  13)  to  35%  (models  71  through  73).  Larger  Pr(l]  yields 
more  frequent  recurrence  of  states  1  and  -1,  which  means  more  elevation  variations  and  hence  a 
rougher  topography.  Intuitively  any  arbitrary  sequence  of  0,  1,  -1  should  suffice  to  generate  a 
corresponding  random  profile.  In  our  LFD  calculations,  however,  the  profile  needs  to  be  sampled  well 
enough  so  that  each  segment  consists  of  at  least  two  subintervais  with  the  same  slope.  As  a  result,  the 
“actual"  probability  of  states  1  and  -1  is  slightly  smaller  than  that  of  the  input.  It  is  interesting  to  note  that 
ft,  happens  to  represent  the  “actual"  probability  of  states  1  and  -1  combined,  which  is  slightly  smaller 
than  twice  that  of  Prfl].  Most  statistics  as  defined  above  show  an  increasing  trend  with  the  increasing 
Pr[l  ],  except  ft4  which  seems  to  exhibit  no  clear  pattern  at  all. 


Table  1.  Roughness  Statistics  of  27  Topographic  Profiles 


Model 

Oi 

«2 

«3 

a* 

fls 

fle 

«7 

Pffl] 

Seed 

11 

0.091 

0.0178 

0.215 

0.0198 

0.240 

0.839 

1.038 

e m 

22887 

12 

0.097 

0.0186 

0.487 

0.0041 

0.511 

0.050 

22892 

13 

0.097 

0.0190 

0.218 

0.0203 

0.260 

0.967 

1.040 

0.050 

85995 

14 

0.146 

0.0225 

0.234 

0.0266 

0.266 

1.389 

0.070 

91333 

15 

0.144 

0.0228 

0.209 

0.0329 

0.310 

1.439 

0.070 

91350 

16 

0.127 

0.0218 

0.265 

0.0181 

0.209 

1.194 

1.053 

0.070 

91285 

21 

0.180 

0.0245 

0.686 

0.0038 

0.528 

1.600 

1.075 

0.100 

93337 

22 

0.182 

0.0248 

0.557 

0.0059 

0.638 

1.757 

1.076 

0.100 

93347 

23 

0.150 

0.0231 

0.435 

0.0079 

0.370 

1.429 

1.062 

0.100 

93356 

31 

0.310 

0.0323 

1.329 

0.0018 

0.981 

2.357 

1.128 

0.150 

748546615 

32 

0.230 

0.0283 

0.234 

0.0421 

0.317 

1.742 

1.095 

0.150 

20333 

33 

0.233 

0.0282 

0.408 

0.0140 

0.303 

2.192 

1.096 

0.150 

86016 

41 

0.283 

0.756 

0.0049 

0.599 

2.357 

1.117 

0.200 

93397 

42 

0.298 

0.0318 

0.878 

0.0039 

0.751 

2.508 

1.124 

0.200 

93407 

43 

0.316 

0.0324 

0.751 

0.0056 

0.819 

2.392 

1.131 

0.200 

12345 

51 

0.401 

0.0366 

0.469 

0.0183 

0.584 

2.609 

1.166 

0.250 

21691 

52 

0.361 

0.0348 

0.553 

0.0118 

0.549 

2.742 

1.150 

0.250 

21696 

53 

0.399 

0.0364 

0.506 

0.0156 

0.633 

2.591 

1.165 

0.250 

86037 

61 

0.478 

0.0400 

1.003 

0.0048 

1.741 

2.676 

1.198 

0.300 

93459 

62 

0.437 

0.0386 

1.160 

0.0032 

1.153 

2.662 

1.181 

0.300 

93468 

63 

0.438 

0.0383 

0.745 

0.0079 

1.112 

2.274 

1.181 

0.300 

93479 

64 

0.477 

0.0400 

0.539 

0.0164 

0.663 

2.685 

1.197 

0.320 

55473 

65 

0.483 

0.0405 

0.934 

0.0055 

1.207 

2.585 

1.200 

0.320 

55499 

66 

0.512 

0.0421 

0.604 

0.0140 

0.795 

2.699 

1.212 

0.320 

55446 

71 

0.502 

0.0410 

0.538 

0.0174 

0.749 

2.475 

1.208 

0.350 

22192 

72 

0.595 

0.0444 

0.728 

0.0112 

0.777 

1.993 

1.246 

0.350 

31087 

73 

0.555 

0.0431 

0.863 

0.0075 

1.129 

2.341 

1.230 

0.350 

86061 

1 
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2.3  Correlation  Between  Various  Roughness  Measures 

The  correlation  coefficients  between  various  roughness  measures  are  listed  in  Table  2.  As  would 
be  expected  from  Table  1 ,  almost  all  of  the  parameters  are  positively  correlated  with  each  other  with  the 
exception  of  Q4.  Thus  it  may  be  anticipated  that,  except  for  Q*  a  profile  which  is  “rough"  in  one  sense 
would  probably  be  rough  in  another  sense  as  well.  The  question  is,  then,  which  parameter  can  optimally 
characterize  the  wave  propagation  phenomenon.  That  is,  which  one  of  the  roughness  parameters  is  the 
best  statistic  to  use  in  predicting  the  amplitude  decay  when  some  previously  uncalibrated  rugged  paths 
are  involved.  The  answer  to  such  a  question  could  well  be  dependent  on  the  seismic  phases  of  interest, 
since  different  seismic  phase  might  interact  with  the  rough  topography  differently.  Throughout  this  study, 
Rg  wave  packet  is  injected  into  the  model  as  the  incident  source,  and  the  transmitted,  reflected  and  the 
scattered  synthetics  are  then  analyzed.  The  ultimate  output  of  this  exercise  would  be  some  calibration 
curves  which  may  be  suitable  for  predicting  the  Rg  amplitudes,  given  an  arbitrary  rugged  profile.  The 
same  procedure  should  be  repeated  for  other  crustal  phases  to  establish  the  associated  relationship  with 
the  roughness  measures,  which  will  be  addressed  in  the  follow-up  report. 


Table  2.Correlation 

Coefficients 

Between  Roue 

ihness  Param 

eters 

«3 

O4 

Os 

O* 

O7 

Prfl] 

Qi 

0.99 

0.58 

-0.27 

0.73 

0.82 

1.00 

0.99 

q2 

_ 

0.60 

-0.27 

0.73 

0.87 

0.99 

0.99 

Oa 

-0.76 

0.84 

0.64 

0.58 

0.55 

_ 

_ 

-0.58 

-0.35 

-0.27 

-0.27 

Os 

_ 

_ 

_ 

_ 

0.67 

0.73 

0.72 

Oe 

0.82 

0.84 

O7 

0.99 
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3.  NUMERICAL  EXPERIMENTS 


3.1  Finite-difference  Scheme 

The  2-0  linear  finite-difference  (LFO)  scheme  popularized  by  Kelly  et  al.  (1976)  is  used  with  the 
Clayton  and  Engquist  (1977)  (Emerman  and  Stephen,  1983)  absorbing  boundary  conditions  on  both  the 
sides  and  the  bottom  of  the  grid.  Even  with  the  absorbing  boundary  conditions,  the  simulations  are  run 
with  large  enough  grids  such  that  for  the  duration  of  the  simulations  reflections  from  the  sides  do  not 
seriously  contaminate  the  results.  The  free-surface  boundary  conditions  follow  from  the  work  of  Jih 
et  al.  (1988)  and  Jih  (1993a)  which  allows  for  the  introduction  of  arbitrary  polygonal  topography  into  the 
wave  propagation  simulations  and  which  has  been  applied  by  McLaughlin  and  Jih  (1988)  for  the  model¬ 
ing  of  teleseismic  P  waves  from  explosive  line  sources  under  realistic  topographic  surfaces.  These 
boundary  conditions  are  extensions  to  those  of  llan  (1975  and  1977).  Recently  a  number  of  techniques 
have  been  pursued  in  an  effort  to  improve  the  computational  performance  of  LFD  solutions  to  wave 
equations.  These  include  variable  grids,  higher  order  schemes,  or  implicit  rather  than  explicit  methods. 
Jih  et  al.  (1988)  have  found  that  none  of  these  methods  provide  significant  improvement  over  the  tradi¬ 
tional  explicit  second  order  centered  LFD  schemes  described  by  Alterman,  Boore,  and  Kelly  when  rough 
topographic  profiles  are  used  at  the  free  surface. 

In  the  present  experiments,  the  medium  has  P-  and  5- wave  velocities  of  5.02  and  2.898  km/sec 
respectively,  corresponding  to  a  Poisson's  ratio  of  0.35  and  a  Rayleigh  wave  velocity  of  2.665  km/sec. 
The  formulation  generating  the  incident  pure  Rayleigh  wave  packet  is  analogous  to  Boore  (1970)  or 
Munasinghe  and  Famell  (1973)  which  has  been  used  extensively  by  Martel  et  al.  (1977),  Fuyuki  and 
Matsumoto  (1980),  Fuyuki  and  Nakano  (1984),  Levander  (1985),  and  McLaughlin  and  Jih  (1986  and 
1987).  The  initial  waveform  is  the  fundamental  mode  Rayleigh  wave  in  the  homogeneous  half  space 
with  a  Ricker  wavelet  shape.  This  method  is  superior  to  another  approach  which  shakes  several  mesh 
points  on  the  free  surface  or  the  grid  edge  to  generate  the  incident  Rayleigh  wave,  since  the  latter  would 
excite  undesirable  body  waves  as  well.  The  pure  Rayleigh  wave  packet  of  Ricker-wavelet  shape  is  par¬ 
ticularly  suitable  for  the  numerical  modeling  study,  since  it  is  well  localized  in  both  the  time  and 
wavenumber  spaces.  To  reduce  the  grid  dispersion,  we  have  oversampled  the  incident  wave  packet  by 
choosing  the  dominant  wavelength  to  be  3.2  km  =  32  AX,  where  AX  =0.1  km  is  the  grid  spacing  in  the 
horizontal  direction.  The  dominant  wavenumber  is  2.358  radians/km  which  corresponds  to  a  wavelength 
of  2.665  km  at  1Hz.  Thus  the  energy  of  our  incident  Rg  wave  packet  is  centered  at  1  Hz. 

The  grids  used  in  LFD  simulations  are  60  km  in  length  and  50  km  in  depth,  with  a  grid  spacing  of 
0.1  km  in  both  X  and  Z  directions.  The  initial  incident  waveform,  the  reflected  and  the  transmitted 
waveforms  were  then  measured  at  free-surface  “sensors".  Transmission  coefficients  and  the  attenua¬ 
tion  factor  Q(f)  were  then  computed  in  the  frequency  domain.  The  useful  bandwidth  is  limited  to  fre¬ 
quencies  below  2.4  Hz  due  to  the  limitation  of  2nd  order  LFD  scheme.  Although  we  quote  dimensions 
in  km  and  velocities  in  km/s,  it  is  possible  to  scale  these  simulations  to  higher  or  lower  frequencies  as 
desired.  As  explained  earlier,  the  LFD  scheme  must  discretize  the  topography  into  well-sampled 
polygons,  which  is  exactly  the  reason  for  adopting  a  Markov  sequence  rather  than  a  white  sequence. 
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Therefore  the  smallest  scale  length  erf  the  topography  is  always  twice  the  horizontal  grid  spacing.  The 
largest  scale  length  is  limited  by  the  total  span  of  the  rugged  profile  used  in  the  simulation. 

3.2  Half-space  Experiment 

Experiments  were  first  conducted  to  justify  that,  in  the  case  of  a  homogeneous  half-space,  the  grid 
dispersion  only  slightly  modulates  the  spectra  of  the  propagating  Rayleigh  wave.  A  particularly  informa¬ 
tive  way  to  visualize  wave  propagation  through  the  models  is  by  using  numerical  Schlieren  diagrams. 
These  “snapshots"  of  the  wavefieid  traveling  through  the  grid  are  generated  by  saving  the  displacement 
(or  strain)  at  each  grid  point  for  a  given  time  step.  Figure  10  shows  the  snapshots  of  the  displacement 
field  at  time  incidents  of  4  seconds  apart.  Shading  in  each  snapshot  is  proportional  to  V  Uz  +  W2  where 
U  and  W  are  the  displacement  along  the  X  and  Z  directions,  respectively,  of  the  grid.  The  darkness  of 
the  plot  is  normalized  to  the  maximum  in  each  frame.  No  significant  attenuation  of  the  Rayleigh  wave 
could  be  measured  for  the  half-space  in  the  frequency  band  of  interest  (Figure  11). 

3.3  Experiments  with  Rough  Topography 

Following  this  check,  we  conducted  27  simulations  with  the  random  topographic  profiles  listed  in 
Table  1.  Each  simulation  took  4000  time  steps.  As  in  the  half-space  case,  Figures  12  through  38  show 
snapshots  of  the  displacement  field  due  to  a  Rayleigh  wave  incident  upon  various  topographic  profiles. 
Successive  frames  are  separated  by  intervals  of  4  seconds,  and  the  incident  Rayleigh  wave  always  trav¬ 
els  from  left  to  right.  In  general,  the  diffraction  pattern  becomes  much  more  complicated  when  the  the 
topography  is  rougher.  This  is  obvious  in  comparing  Figure  13  (for  model  No.12)  and  Figure  32  (for 
model  No.63),  for  instance.  In  all  cases  the  comers  act  as  point  sources  radiating  converted  body  waves 
(mainly  S  wave). 

Thirty  sensors  numbered  from  No.1  through  No.30  are  placed  2  km  apart  along  the  free  surface  to 
record  the  displacement.  Figures  39  through  66  show  the  seismic  section  of  the  vertical  component  of 
each  model  with  the  maximum  zero-to-peak  amplitude  shown  on  the  right.  The  distance  between  each 
sensor  and  the  “epicenter”  of  the  initial  Rayleigh  wave  packet  is  also  shown.  Smoother  models  tend  to 
have  longer  reflected  wave  trains  than  do  rougher  models.  Also,  the  peak  amplitude  of  the  transmitted 
signal,  which  is  measured  with  trace  No.29  of  each  model,  drops  significantly  as  the  model  becomes 
rougher.  The  reflected  waves  systematically  consist  of  multiply  reflected  Rayleigh  waves  which  last  22 
seconds  or  longer,  which  is  the  time  required  for  the  Rg  packet  to  cross  the  whole  profile  and  then 
bounce  back.  On  the  other  hand,  the  transmitted  synthetic  is  dominated  by  very  few  cycles  of  Rayleigh 
waves.  In  the  following,  these  qualitative  observations  are  further  elaborated  in  a  quantitative  manner. 

The  power  spectra  of  traces  No.2  and  No.29  of  each  model  are  measured  to  represent  the 
reflected  and  transmitted  signal.  We  define  the  Rayleigh  wave  energy  transmission  coefficient,  T,  as 
T(f)  =  Power, rans(f)/Powerjnc(f)  where  Power^f)  is  the  power  spectrum  of  the  incident  Rayleigh  wave. 
The  reflection  coefficient,  R,  is  defined  in  an  analogous  way:  R(f)  =  Powerr9fl(f)/Powerinc(f).  Table  2 
gives  the  statistics  of  all  27  topographic  profiles  used  and  the  corresponding  attenuation  factor  Q  at  1 
Hz,  which  is  the  dominant  frequency  of  the  incident  Rg  packet.  Here  the  quality  factor,  Q,  is  computed 
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via  y  ■  Q_1(f)  ■  -ln(T(f))  -  (V^idX)  as  in  McLaughlin  and  Jih  (1986  and  1987). 


Table  3.  TransmisstorVRefl< 

act  ion/ Sc 

altering  Sts 

atlstics 

Model 

T(1Hz) 

R(1Hz) 

1-T-R 

Q(1  Hz) 

Y 

11 

0.360 

0.181 

0.459 

65.6 

0.0153 

12 

0.345 

0.251 

0.404 

72.0 

0.0139 

13 

0.371 

0.207 

0.422 

81.1 

0.0123 

14 

0.174 

0.198 

0.628 

40.4 

0.0248 

15 

0.320 

0.207 

0.473 

58.2 

0.0172 

16 

0.284 

0.225 

0.491 

52.4 

0.0191 

21 

0.138 

0.237 

0.625 

33.9 

0.0295 

22 

0.217 

0.289 

0.494 

50.5 

0.0198 

23 

0.148 

0.188 

0.664 

37.1 

0.0269 

31 

0.041 

0.130 

0.829 

22.3 

0.0449 

32 

0.076 

0.249 

0.675 

28.2 

0.0354 

33 

0.189 

0.193 

0.618 

46.4 

0.0215 

41 

0.063 

0.181 

0.756 

26.0 

0.0385 

42 

0.049 

0.160 

0.791 

27.0 

0.0371 

43 

0.165 

0.143 

0.692 

44.0 

0.0227 

51 

0.084 

0.101 

0.815 

31.3 

0.0320 

52 

0.047 

0.229 

0.724 

0.0439 

53 

0.189 

0.256 

0.555 

44.3 

0.0226 

61 

0.209 

0.073 

0.718 

64.9 

0.0154 

62 

0.123 

0.112 

0.765 

37.6 

0.0266 

63 

0.053 

0.064 

0.883 

23.7 

0.0422 

64 

0.051 

0.899 

22.2 

0.0450 

65 

0.133 

0.079 

0.788 

35.1 

0.0285 

66 

0.066 

0.043 

0.891 

26.3 

0.0380 

71 

0.091 

0.068 

0.841 

32.3 

0.0309 

72 

0.011 

0.059 

0.930 

16.8 

0.0596 

73 

0.006 

0.047 

0.947 

15.2 

0.0658 

A  significant  conclusion  is  immediately  apparent  from  Table  3,  that  the  reflected  energy  and 
transmitted  energy  is  far  from  100%  of  the  incident  energy,  and  that  body  wave  conversion  must  be  very 
efficient.  Even  the  smoothest  model  (such  as  model  No.  12)  caused  a  scattering  loss  of  at  least  40%.  In 
the  extreme  case  of  model  No.73,  nearly  95%  of  the  Rayleigh  wave  energy  is  lost  to  body  wave  conver¬ 
sion,  which  results  in  a  very  low  quality  factor.  The  model  which  yields  the  highest  value  of  Q  has  a 


transmission  coefficient  of  about  37%,  and  about  42%  of  Rg  energy  is  lost  to  scattering.  The  spatial  Q 
values  associated  with  rough  topography  as  listed  in  Table  3  are  in  the  same  range  as  the  Op  of  the 
uppermost  crust  in  many  places  including  the  Arabian  shield  (Mokhtar,  et  a/.,  1988),  the  eastern  folded 
region  of  the  Arabian  Peninsula  (Seber  and  Mitchell,  1991),  Maine  (Al-Eqabi  and  Herrmann,  1993),  and 
different  portions  of  New  England  (Safcia  et  at.,  1990;  Saikia,  1992). 

The  nature  of  the  reflected  and  transmitted  signals  can  be  determined  with  the  polarity.  Figures  68 
through  94  show  the  synthetics  at  surface  sensors  No.2  and  No.29,  which  have  been  used  to  represent 
the  reflected  and  transmitted  waves,  respectively,  in  Table  3.  Vertical  and  horizontal  components  are 
plotted  in  solid  and  dashed  lines,  respectively.  The  particle  motion  shown  on  the  left  dearly  exhibits  the 
retrograde  motion.  For  comparison,  the  synthetics  recorded  at  a  deep  sensor  are  shown  at  the  bottom 
of  each  figure.  This  sensor,  being  37.5  km  in  depth,  records  only  the  scattered  body  waves,  as  is  indi¬ 
cated  by  the  nearly  linear  partide  motions. 

The  four  propagation  parameters,  namely  y.  T,  R,  and  1-T-R,  have  been  regressed  against  each 
individual  roughness  coefficients  separately,  as  shown  in  Figures  95  through  103.  Table  4  summarizes 
the  correlation  between  various  roughness  measures  and  the  Rg  propagation  characteristics,  £2,  and 
£22  have  the  highest  correlation  (0.87  and  0.88,  respectively)  with  the  scattering  loss.  This  is  also  illus¬ 
trated  in  lower-right  plot  of  Figures  95  through  103.  Note  that  £2,  and  Q2  perform  better  in  predicting  the 
scattering  loss  than  in  predicting  the  transmission  or  reflection  alone.  For  a  practical  monitoring  situation, 
however,  it  might  be  desirable  to  infer  the  transmission  coefficient,  T,  of  an  uncatibrated  path  in  order  to 
gain  some  information  about  how  the  seismic  wave  behaves.  On  this  regard,  £21t  £2Z,  and  £2*  might 
serve  the  purpose  equally  well,  with  the  correlation  coefficients  around  0.73-0.76.  04,  the  Durbin- 
Watson  number  or  the  von  Neumann  ratio,  turns  out  to  be  the  poorest  statistics  in  all  tests.  The  correla¬ 
tion  coefficient  between  £24  and  any  of  the  propagation  statistics  is  always  no  higher  than  0.27  (Figure 
98).  As  discussed  earlier,  £24  is  essentiaty  the  square  of  the  quotient  of  two  other  roughness  measures, 
namely  £22  and  £23.  Thus  combining  these  two  parameters  offers  no  advantage  at  all.  That  £2«  is  nega¬ 
tively  correlated  to  all  other  roughness  measures  {cf.  Table  2),  scatter  loss,  as  well  as  y,  would  suggest 
that  £24  is  not  a  good  candidate  for  quantifying  the  roughness.  Several  systematic  patterns  can  be 
observed  from  the  regression  results;  both  y  and  the  scattering  loss  (1-T-R)  increase  with  the  roughness 
coefficients  (except  £24).  whereas  T  and  R  decrease  with  the  roughness  coefficients. 


Table  4.Correlation  Between  Roughness  and  Propagation 

Characte 

iristics 

T(1Hz) 

R(1Hz) 

1-T-R 

£J2 

^3 

«5 

«6 

a7 

P'11] 

Y 

-0.88 

-0.54 

0.84 

0.67 

0.66 

0.42 

-0.17 

0.31 

0.49 

0.66 

0.64 

T(1Hz) 

0.54 

-0.92 

-0.73 

-0.76 

-0.54 

0.24 

-0.43 

-0.74 

-0.73 

-0.73 

R(1Hz) 

-0.83 

-0.81 

-0.79 

-0.49 

0.24 

-0.65 

-0.57 

-0.81 

-0.80 

1-T-R 

...  .  _ 

0.87 

0.88 

0.59 

•0.27 

0.59 

0.76 

0.86 

0.86 
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Table  5.  Regression  Results 


Slope 


0.057±0.013 


-0.501  ±0. 


-0.39310. 


0.89410.104 


1 .06910.242 


-9.91511.686 


-7.20611.117 


17.12111.869 


0.01910. 


-0.19910. 


-0.12710.045 


0.32610. 


Intercept 


0.01310.004 


0.30610. 


<*Y) 


0.010 


0.077 


T(1Hz) 


R(1Hz) 


1-T-R 


0.41310. 


-0.00310. 


0.46210. 


0.38410.037 


0.15410.061 


0.01910. 


.041 


T(1Hz) 


R(1Hz) 


1-T-R 


0.010 


0.072 


T(1Hz) 


R(1Hz) 


1-T-R 


0.013 


0.094 


0.135 


04 

1 

•0.23410.274 

0.03410.004 

0.014 

-.17 

04 

T(1Hz) 

2.66712.181 

0.24 

O4 

R(1Hz) 

1.87811.527 

0.13210.024 

0.076 

0.24 

«4 

1-T-R 

-4.54613.239 

0.75410.052 

0.161 

-.27 

T(1Hz) 


R(1Hz) 


1-T-R 


0.01110.007 


-0.12710.054 


-0.13610.032 


0.26310.072 


0.01110.004 


-0.13110.024 


-0.07010.020 


0.20210.035 


.041 


.24810.024 


0.51810. 


0.013 

0.31 

0.101 

-.43 

0.059 

-.65 

0.135 

0.59 

0.012 

0.49 

0.075 

-.74 

0.064 

-.57 

0.109 

0.76 

T(1Hz) 


R(1Hz) 


1-T-R 


0.41810.051 


0.30110.044 


0.28110.074 


O7 

1 

0.13710.031 

-0.12410.035 

0.010 

0.66 

O7 

T(1Hz) 

-1.210+0.229 

1.51710.259 

0.077 

-.73 

o7 

R(1Hz) 

-0.94910.136 

1.22910.154 

0.046 

-.81 

o7 

1-T-R 

2.15910.251 

-1.74610.284 

0.084 

0.86 

Prtl] 

y 

0.08010.019 

0.01510.004 

0.011 

0.64 

Prfl] 

T(1Hz) 

-0.73410.139 

0.29410.031 

0.077 

-.73 

Prtl] 

R(1Hz) 

-0.56610.085 

0.26910.019 

0.047 

-.80 

Prtl] 

1-T-R 

1.300+0.156 

0.43710.035 

0.086 

0.86 

3.5  Scattered  Body  Waves 

A  horizontal  array  is  formed  of  30  sensors  spaced  2  km  apart  at  a  depth  of  37.5  km  to  record  the 
converted  P  and  SV  wavefieids.  Figures  103  through  129  give  the  seismic  sections  of  the  converted  P 
wave  (top)  and  S  wave  (bottom)  displacements  for  each  model.  The  take-off  angle  as  well  as  the  zero- 
to-peak  amplitude  of  each  syntactic  are  shown  on  the  right.  The  seismic  sections  are  plotted  in  a  uni¬ 
form  scale.  Note  that  the  converted  SV  waves  typicaly  have  a  peak  amplitude  two  to  three  times  larger 
than  that  of  the  the  converted  P  waves.  The  P  and  S  wavefieids  can  be  seen  to  be  radiated  from 
numerous  point  diffractors  along  the  topographic  profile,  as  was  already  suggested  by  the  snapshots 
(Figures  12-38).  For  the  lower  portion  of  the  seismic  section,  i.e. ,  synthetics  numbered  50  and  on,  the 
take-off  angles  are  generally  more  oblique  than  30°.  If  a  velocity  gradient  were  introduced  into  the 
model,  most  of  these  rays  would  pass  the  critical  angle  and  could  well  be  multiply  reflected  in  the  crustal 
waveguide. 

In  plotting  the  seismic  sections,  each  synthetic  seismogram  is  steered  towards  its  ray  direction 
which  is  determined  with  the  particle  motion  within  an  appropriate  time  window.  To  illustrate  the  robust¬ 
ness  of  this  approach  in  separating  P  and  SV  components,  Figures  130  through  156  display  three  syn¬ 
thetic  seismograms  for  each  model.  The  particle  motion  throughout  the  whole  LFD  simulation  (i.e. ,  40- 
second  window)  is  shown  on  the  left,  with  the  P  and  SV  components  shown  in  solid  and  dashed  lines, 
respectively,  on  the  right.  Traces  31,  40,  and  59  of  each  model  are  plotted  from  top  to  bottom.  The 
take-off  angle  of  each  ray  is  determined  by  the  P  onset,  and  it  turns  out  that  within  the  initial  P  window 
(typically  from  10  to  15  seconds),  the  SV  component  is  nearly  negligible.  On  the  other  hand,  the  later 
arrivals  (e.g. ,  20  seconds  and  on)  are  dominated  by  SV  component  with  peak  amplitude  a  few  times 
larger  than  that  of  the  P  waves,  except  for  the  obliquely  scattered  trace  (e.g.,  the  synthetic  No.59)  of  a 
few  models. 

Although  the  topography  is  embedded  in  the  central  portion  of  the  grid,  the  scattered  5  wave  field 
is  far  from  isotropic.  The  backscattered  S  waves  have  steeper  take-off  angles  and  longer  duration, 
whereas  the  forward  scattered  S  waves  have  more  compact  wavetrains.  There  are  twice  as  many  sen¬ 
sors  that  recorded  forward-scattered  S  waves.  Furthermore,  many  of  them  have  very  oblique  take-off 
angles.  This  is  in  support  of  the  interpretation  that  we  are  observing  a  series  of  S-wave  point  sources 
spread  along  the  free  surface  which  are  “triggered”  in  the  forward  direction.  Thus  even  if  we  assume 
that  the  subsurface  geology  is  perfectly  homogeneous  and  isotropic,  the  scattered  S  wavefield  will  not 
be  as  symmetric.  Both  the  forward  and  backward  directions  have  about  the  same  number  of  rays  that 
would  penetrate  into  the  deeper  earth.  However,  there  is  always  a  significant  amount  of  energy  which  is 
favorably  trapped  in  the  crustal  waveguide  in  the  forward  direction. 
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4.  DISCUSSION  AND  CONCLUSIONS 


In  this  study,  we  have  performed  2-D  numerical  Rayleigh-wave  scattering  experiments  upon  com¬ 
plicated  random  topographies.  The  simulations  show  that  reflection  is  an  efficient  process  and  that  the 
bufc  of  the  energy  which  is  not  transmitted  across  the  topography  is  either  converted  to  body  waves  or 
reflected  back.  Through  the  simulations,  we  illustrate  that  P  and  S  waves  are  radiated  as  coda  to 
teleseismic  distances,  and  some  P-SV  energy  is  trapped  in  the  crustal  waveguide  as  P0  ,  Sg ,  and  Lg . 
The  body  waves  appear  to  radiate  from  the  changes  in  slope  and  to  first  order  can  be  visualized  as 
being  from  point-force  diffractors.  Mi  these  observations  are  consistent  with  our  previous  study 
(McLaughlin  and  Jih,  1986  and  1987;  Jih  and  McLaughlin,  1988). 

The  spatial  Q  values  associated  with  rough  topography  obtained  in  the  study  are  in  the  same 
range  as  the  Op  of  the  uppermost  crust  in  many  places  including  the  Arabian  shield  (Mokhtar,  et  at. , 
1988),  the  eastern  folded  region  of  Arabian  Peninsula  (Saber  and  Mitchell,  1991),  and  Maine  (AJ-Eqabi 
and  Herrmann,  1993).  Saikia  (1992)  (and  Saikia  et  at. ,  1990)  also  finds  that  a  low  Op  of  15  is  needed 
in  order  to  match  his  Rg  synthetics  with  the  data  observed  at  southwestern  New  England.  An  accurate 
characterization/prediction  of  the  propagation  of  regional  phases  due  to  explosion  sources  in  regions  of 
high  proliferation  concern  requires  an  accurate  regionalized  map  of  the  attenuation  coefficient.  Although 
a  global  Lg  coda  Q  map  is  available  for  much  of  Eurasia  and  some  other  regions  of  the  world,  the 
identification  of  small  seismic  events  has  raised  a  need  to  refine  such  a  map  for  a  more  detailed  region¬ 
alization  for  complex  areas.  It  is  known  that  changes  in  the  quality  factor  of  the  propagation  path  can 
have  a  significant  influence  on  magnitude  estimates  and  therefore  on  network  performance  estimates 
(e.g. ,  Jih  and  Lynnes,  1993).  Regionalized  Q  maps  (for  different  phases)  need  to  be  established  and 
must  be  coupled  with  the  appropriate  source  model  tor  the  distinct  geologies. 

A  major  issue  for  the  Non-Proliferation  Treaty  is  the  discrimination  of  large  chemical  explosions 
from  possible  clandestine  or  small  nuclear  tests.  Unless  discrimination  is  possible,  the  numerous  mining 
blasts  could  give  ample  opportunity  for  conceal ng  clandestine  tests.  Rippled-fired  explosions  are  com¬ 
monly  used  to  fragment  rocks  during  quarry  and  open-pit  mining.  Through  LFD  calculations,  Jih  (1993b) 
finds  that  ripple  firing  could  excite  directionally  dependent  Rg  and  S*.  Thus  the  lack  of  Rg  may  not  be 
always  indicative  of  a  deep  source.  Rather,  it  could  also  be  due  to  the  shot  pattern.  However,  the 
enhanced  Rg  in  the  forward  direction  of  ripple  firing  can  be  strongly  attenuated  by  lateral  heterogeneity 
and  surface  topography,  as  is  illustrated  in  the  present  study.  The  scattered  Rg  energy  could  then  cou¬ 
ple  into  the  crustal  waveguide  as  Lg  and  other  phases.  Since  such  scattering  mechanisms  are  com¬ 
monly  present  in  many  quarry  sites  or  mines,  it  is  not  surprising  that  the  directional  enhancement  of  Rg 
may  not  be  always  observable.  The  spall  could  obscure  the  azimuthal  dependency  of  Rg  as  well.  LFD 
modeling  studies  suggest  that  the  Rg  -to -S  scattering  is  in  general  stronger  than  that  of  Rg  -to -P,  which 
could  provide  a  plausible  explanation  of  why  quarry  blasts  and  mining  blasts  should  discriminate  less 
well  from  earthquakes  than  would  contained  nuclear  explosions. 

The  Lg  phase  has  an  increasingly  prominent  role  in  nuclear  test  detection,  discrimination,  and 
yield  estimation.  Regional  Lg  phase  is  a  short-period  guided  wave,  composed  mainly  of  a  sequence  of 
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multiply  reflected  posterities  S- waves  trapped  in  the  crustal  waveguide.  The  formation  of  Lg  requires  a 
significant  amount  of  energy  propagating  as  S  waves,  either  directly  radiated  from  the  source  or  con¬ 
verted  from  other  phases,  and  these  S  waves  must  be  properly  trapped  in  the  crustal  wave  guide  (Xie 
and  Lay,  1993;  Frankel,  1989;  and  others).  Shallow  explosions  are  strong  Rg  generators,  and  these 
explosion-generated  Rg  waves  can  contribute  to  Lg  via  scattering  at  rough  topohraphy.  LFD  simulations 
using  a  pure  Rg  packet  as  the  initial  condition  offers  a  promising  approach  to  quantify  effects  of  scatter¬ 
ing  by  topography  and  heterogeneity  upon  surface  waves,  as  illustrated  in  this  study.  A  complete 
exploration  of  the  problem  will  require  variation  of  the  P  and  S  velocities,  near  surface  velocity  gra¬ 
dients,  as  well  as  crustal  velocity  heterogeneity,  since  the  bw-Q  sediments  can  be  equally  attenuative, 
as  reported  by  Mitchell  and  Hwang  (1987).  Adding  any  of  these  mechanisms  to  the  numerical  experi¬ 
ments  can  only  increase  the  S  coupling  ( cf .  Aki,  1992).  Thus  we  can  anticipate  that  in  a  more  realistic 
geological  environment,  the  Rayleigh-to-L0  conversion  would  be  more  profound  than  in  the  case  of  a 
half  space.  Recently  Patton  and  Taylor  (1993)  observed  Lg  spectral  null  around  0.55  Hz  for  Yucca  Flat 
explosions,  which  could  be  due  to  the  null  in  excitation  of  Rg  by  a  CLVD  source  as  a  function  of  depth. 
The  null  can  not  be  accounted  for  by  the  delay  between  P  and  spall  slapdown  when  those  differ  for 
shots  at  constant  depth  but  different  yield.  Thus  the  scattering  of  large  amplitude  Rg  waves  in  the  near 
field  may  be  the  most  important  source  of  Lg  waves  with  frequencies  near  1  Hz  for  explosions  on  Yucca 
Flat  or  mountainous  site. 

Zhang  and  Lay  (1993)  found  the  Sn/Lg  amplitude  ratio  to  be  strongly  dependent  on  some  rough¬ 
ness  measures  of  paths  from  Semipalatinsk  to  several  Eurasian  stations.  Perhaps  S„  interacts  with  the 
irregular  Moho  topography  in  a  way  very  similar  to  that  of  Rg  with  the  free-surface  topography.  Thus 
this  may  suggest  that  energy  partitioning  in  the  regional  wavefield  is  controlled  by  waveguide  structure 
and  attenuation  variations  which  are  manifested  in  the  surface  topography.  Further  investigation  along 
this  line  with  numerical  modeling  techniques  as  presented  in  this  study  would  be  very  useful.  A  distinct 
advantage  of  using  a  source-free  Rayleigh  wave  packet  in  the  simulation  (as  adopted  in  this  study)  is 
that  the  interaction  with  certain  types  of  heterogeneity  (such  as  the  topography)  can  be  isolated.  Scaling 
our  results  with  other  source-dependent  relationships  (such  as  the  depth-dependence  of  Rg  excitation)  is 
straightforward. 

The  simple  functional  relationship  between  some  of  the  roughness  parameters  might  lead  us  to 
regard  some  of  the  roughness  parameters  as  redundant.  However,  it  turns  out  that  the  performance  of 
these  roughness  measures  could  be  different.  This  simply  means  that  by  applying  appropriate  “transfor¬ 
mation”  on  the  data,  we  might  be  able  to  find  a  better  linear  regression.  A  useful  family  of  transforma¬ 
tions  for  this  purpose  is  the  family  of  Box-Cox  transformations  (see,  for  example,  Johnson  and  Wichern, 
1992).  When  we  regress  y  and  the  transmission  power  (T(1Hz))  on  any  of  the  roughness  measures, 
say  Q1(  it  should  not  be  surprising  to  see  that  a  linear  fit  is  appropriate  for  y,  but  not  T(1Hz).  Converting 
T(1Hz)  to  y  with  a  simple  logarithmic  operation  is  a  typical  example  of  the  Box-Cox  transformation.  Such 
an  exercise  of  transforming  nonnormal  data  into  more  nearly  “normal  looking”  data  (around  the  calibra¬ 
tion  curve)  is  a  standard  practice  in  the  statistics  field,  and  it  may  be  worth  exploring. 
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7.  FIGURE  CAPTIONS 


Figure  1:  3  rugged  profiles  corresponding  to  models  No.  11  through  No.  13,  respectively,  are  shown  on 
the  left,  with  the  log-amplitude  wavenumber  spectra  of  the  topography  in  the  center,  and  statis¬ 
tics  characterizing  the  topography  listed  on  the  right.  The  log-amplitude  wavenumber  spectra  are 
normalized  by  the  peak  value  of  the  spectra  (which  is  shown  at  the  right  as  K(0)).  The  seed 
number  and  the  preassigned  probability  of  state  1,  which  determine  the  random  topography,  are 
printed  under  the  profile. 

Figures  2-9:  Same  as  Figure  1  except  for  models  No.  14  through  No.73. 

Figure  10  The  snapshots  of  the  displacement  field  at  time  incidents  of  4  seconds  apart.  Shading  in  each 
snapshot  is  proportional  to  V  U2  +  Wz  where  U  and  W  are  the  displacement  along  the  X  and  Z 
directions,  respectively,  of  the  grid.  The  darkness  of  the  plot  is  normalized  to  the  maximum  in 
each  frame. 

Figure  1 1 :  Spectral  ratio  of  the  transmitted  Rayleigh  waves  over  the  incident  Rayleigh  waves  in  a  homo¬ 
geneous  half-space.  Grid  dispersion  only  modulates  the  spectra  slightly.  The  frequency  band 
of  interest  in  this  study  is  around  1  Hz.  No  significant  attenuation  of  the  Rayleigh  wave  could  be 
measured  for  the  half-space  in  the  frequency  band  of  interest. 

Figure  12:  The  snapshots  of  the  displacement  field  due  to  a  Rayleigh  wave  incident  upon  topographic 
profile  No.1 1 . 

Figures  13-38:  Same  as  Figures  10  and  12  except  for  models  No.12  through  No.73.  In  general,  the 
diffraction  pattern  becomes  much  more  complicated  when  the  the  topography  is  rougher.  This 
is  obvious  in  comparing  Figure  13  (for  model  No.12)  and  Figure  32  (for  model  No.63),  for 
instance.  In  all  cases  the  comers  act  as  point  sources  radiating  converted  body  waves  (mainly 
S  wave). 

Figure  39:  The  seismic  section  of  the  vertical  component  of  the  half-space  model.  The  sensors  are 
evenly  distributed  (2  km  apart)  along  the  free-surface.  The  maximum  zero-to-peak  amplitude  of 
each  trace  is  shown  on  the  right. 

Figures  40-66:  Same  as  Figure  39  except  for  models  No.12  through  No.73.  Smoother  models  tend  to 
have  longer  reflected  wave  trains  than  do  rougher  models.  Also,  the  peak  amplitude  of  the 
transmitted  signal,  which  is  measured  with  trace  No.29  of  each  model,  drops  significantly  as  the 
model  becomes  rougher.  In  the  following,  these  qualitative  observations  are  further  elaborated 
in  a  quantitative  manner. 

Figure  67:  The  synthetics  for  the  half-space  model  at  three  different  distances  are  shown  on  the  right 
with  the  vertical  and  horizontal  components  are  plotted  in  solid  and  dashed  lines,  respectively. 
The  particle  motion  is  plotted  at  the  left.  Note  the  clear  retrograde  motion. 

Figure  68:  The  synthetics  at  surface  sensors  No. 2  and  No.29,  which  have  been  used  to  represent  the 
reflected  and  transmitted  waves,  respectively,  in  Table  3.  Same  as  in  Figure  67,  the  vertical 
and  horizontal  components  are  plotted  in  solid  and  dashed  lines,  respectively.  The  reflected 
waves  (top)  is  composed  of  multiply  reflected  Rayleigh  waves.  On  the  other  hand,  the  transmit¬ 
ted  synthetic  (middle)  is  dominated  by  very  few  cycles  of  Rayleigh  waves.  The  particle  motion 
shown  on  the  left  clearly  exhibits  the  retrograde  motion.  For  comparison,  the  synthetics 
recorded  at  a  deep  sensor  are  shown  at  the  bottom  of  each  figure.  This  sensor,  being  37.5  km 
in  depth,  records  only  the  scattered  body  waves,  as  is  indicated  by  the  nearly  linear  particle 
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motions. 


Figures  69-94:  Same  as  Figure  68  except  for  models  No.12  through  No.73. 

Figure  95:  The  four  propagation  parameters  (y.  T,  R,  and  1-T-R)  are  regressed  against  ft,.  The  corre¬ 
lation  of  0.87  between  ft,  and  the  scattering  toss  is  statistically  significant. 

Figures  96-102:  Same  as  Figure  95  except  for  different  roughness  measures,  ft,  and  ft2  perform 
better  in  predicting  the  scattering  loss  than  in  predicting  the  transmission  or  reflection  atone.  ft4, 
the  Durbin-Watson  number  or  the  von  Neumann  ratio,  turns  out  to  be  the  poorest  statistics  in  all 
tests.  The  correlation  coefficient  between  ft4  and  any  of  the  propagation  statistics  is  always  no 
higher  than  0.27  (Figure  98).  Several  systematic  patterns  can  be  observed  from  the  regression 
results:  both  y  and  the  scattering  toss  (1-T-R)  increase  with  the  roughness  coefficients  (except 
ft4),  whereas  T  and  R  decrease  with  the  roughness  coefficients. 

Figure  103:  The  seismic  sections  of  the  converted  P  wave  (top)  and  S  wave  (bottom)  displacements  of 
model  No.1 1 .  The  take-off  angle  and  the  zero-to-peak  amplitude  of  each  synthetic  are  shown 
on  the  right. 

Figures  104-129:  Same  as  Figure  103  except  for  models  No.12  through  No.73.  Note  that  the  converted 
SV  waves  typically  have  a  peak  amplitude  two  to  three  times  larger  than  that  of  the  the  con¬ 
verted  P  waves.  For  the  tower  portion  of  the  seismic  section,  i.e. ,  synthetics  numbered  50  and 
on,  the  take-off  angles  are  generally  more  oblique  than  30°. 

Figure  130:  Three  (rotated)  synthetic  seismograms,  (traces  No.3i  (top),  No.40  (middle),  and  No.59  (bot¬ 
tom)),  generated  with  model  No.11.  The  P  and  SV  components  are  shown  in  solid  and  dashed 
lines,  respectively.  The  particle  motion  throughout  the  whole  40-second  window  is  shown  on  the 
left.  The  take-off  angle  of  each  ray  is  determined  by  the  P  onset. 

Figures  131-156:  Same  as  Figure  130  except  for  models  No.12  through  No.73.  Note  that  within  the  ini¬ 
tial  P  window  (typically  from  10  to  15  seconds),  the  SV  component  is  nearly  negligible.  Later 
arrivals  (e.g. ,  20  seconds  and  on)  are  dominated  by  SV  component  with  the  peak  amplitude  a 
few  times  larger  than  that  of  the  P  waves,  except  for  the  obliquely  scattered  trace  (e.g. ,  the 
synthetic  No.59)  of  a  few  models  . 
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Rg  (1  Hz:2.665km);  Grid  600*500  [dx«0.1km,dt=10ms] 
Profile  00  (seed  00,  Pr(l)-OO) 


Figure  10 


RATIO  vers  1.5 

Wed  Oct  27  10:16:11  1993 

Model  00 

NOISE  POWER  NOT  SUBTRACTED  FROM  SIGNAL  POWER 
S/N  POWER  THRESHOLD  -  2.0  N 
NO  SMOOTHING 


O 


FREQ  BAND  SLOPE  +/- 

Y INTERC 

AV  INTEG  AV  RATIO 

SD 

0.75-1.25 

5.869O-03 

1.1 180-03 

9. 7760-01 

-3.947O-03 

-3.798O-03 

1.040O-03 

0.50-1.50 

1 .0230-03 

I.OIIe-03 

9.907O-01 

-3 .0470-03 

-3.058O-03 

1 .2010-03 

0.01-1.00 

-3.351O-02 

1.422O-02 

1  047a+00 

2.943O-03 

2. 6850-03 

1. 7470-02 

1.00-2.00 

1 .0370-02 

3.236e-03 

9.670O-01 

9.790O-04 

9.942O-04 

4.884O-03 

0.50-3.00 

5.904O-03 

1 ,235e-02 

9.772O-01 

3.3710-04 

1 .6700-04 

5.272O-02 

Figure  11 
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Rg  (1Hz:2.665km);  Grid  600*500  [dx=0.1km,dt*10ms] 
Profile  11  (seed  22887,  Pr(1)=0.05) 


SORT  (V**2+H"2)  000S6C,  1  600  1  200 

X  Tic  Mark  -  1  Uni,  Z Tic  Mark-  t  Unit 


0  000006 +00  0  364226+00 


FD  sUp  400 

SORT(V**2+H**2)  4  000SEC,  1  600  1  200 

XTicMark-  1  Unit;  ZTeMafk-  1  Unit 
0000006+00  0  297606+00 


FO  sup  600 

SORT(V*"2+H**2)  8.000S6C,  1  600  1  200 

XTioMvk-  1  Unit;  Z Tic  Mark  -  1  Unit 


0000006+00  0  296966+00 


FDstsp  1200 

9QRT(V2+W**2)  12.00GS6C,  1  600  1  200 

X  Tic  Mark-  1  Unit;  Z  Tic  Mark  -  1  Unit 


0.000006+00  0.306676+00 


FDstsp  1600 

SORT(V**2+H**Z)  16.000S6C,  1  600  1  200 


X  Tic  Mark-  1  Unit,  Z  Tic  Mark  -  1  Unit 


0.000006+00  0.224476+00 


I  FDstsp  2000 

SORT(V**2+H**2)  20  000S6C.  1  600  1  200 

X  Tic  Mark  -  1  Unit;  2  Tic  Mark  -  t  Unit 


O.OOOOOE+OO  0  205106+00 


Figure  12 
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Rg  (1Hz:2.665km);  Grid  600*500  [dx- 0.1  km,dt=  10ms] 
Profile  12  (seed  22892,  Pr(1)=0.05) 


SORT (V*2+W“*2)  0  000SEC.  1  MO  1  200 

XTieMnk-  1  Unit;  Z  Tie  Marti.  1  Uni 


0.000006+00  0  316066+00 


FDstap  400 

SQRT(V~2+H~2)  4  000SEC,  1  600  1  200 

X  Tic  Mark  -  1  Uni;  Z  Tic  Marti  -  1  Unit 


0.000006+00  0.342166+00 


FDstap  600 

SORT (V“2+H**2)  6.000SEC.  1  600  1  200 

X  Tie  Marti .  1  Uni;  ZTeUark-  1  Uni 


0000006+00  0  351476+00 


FOstap  1200 

SORT (V~2+H**2)  12000S6C.  1  600  1  200 

X  Tie  Mark-  1  Unit;  Z  Tie  Mark  -  1  Unt 


0  000006+00  0.446526+00 


FDstsp  1600 

SORT(V~2+H**2)  16.000S6C,  1  600  1  200 

X Tic  Mark  -  1  Unit;  Z Tic  Mark-  1  Unit 
0.000006+00  0  442046+00 


FOstap  2000 

SORT  (V**2+H**2)  20  000S6C,  1  600  1  200 

X  Tic  Marti  -  1  Unit;  Z  Tic  Mack  -  1  Unit 


0  000006+00  0.434676+00 


Figure  13 


Ft)  Map  0 

SO HT(V-2*H"2>  OOOS6C  1  too  1  aoo 
X  Tc  Marti  -  1  Uirt;  ZTioMaik.  1  Uni 
0  000006*00  0  470906*00 


FDatap  400 

SQRT(V~2*H”2)  4.000S6C.  1  900  1  200 

X  Tie  Marti-  1  Uni;  ZTieMartt-  1  Una 
0.000006*00  0  399996*00 


FOMap  900 

SQFtr(V-2*H”2)  6.000SEC.  1  000  1  200 

X  To  Marti-  1  Unt:  Z  Tic  Marti-  1  Uni 
0.000006*00  0.309716*00 


FDMp  1900 

8QflT(V“2*K"2)  12-OOOSEC,  1  900  1  200 
X  Tie  Marti*  1  Unt;  ZTieMaik-  1  Unit 
0.000006*00  0301 246*00 


FDaaap  1900 

SOHT(V2*H“2)  10.000S6C.  1  600  1  200 

X  Tie  Marti*  1  Unt;  Z  Tie  Marti.  1  Unit 
0.000006*00  0  295216*00 


FDalap  2000 

SOFtT(V~2*H”2)  20.000SEC.  1  600  1  200 

X  Tic  Marti  .  1  Uni;  Z  Tic  Marti.  1  Unit 
0000006*00  0.287706*00 


Rg  (1Hz:2.665km);  Grid  600*500  [dx=0.1km,dt=10ms] 
Profile  13  (seed  85995,  Pr(1)=0.05) 


Figure  14 
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FDclap  0 

SQRT(V 2+H**2)  00086C,  1  600  1  200 

XTic  Uttfc  •  1  Uni;  Z  Tm  MaA  .  1  Uni 
0000006+00  0  551236+00 


FOtlap  400 

SQHT(V“2+W2)  4.000SEC.  1  000  1  200 

XTfeMaA-  1  Unit;  ZTeMaA-  1  Ur* 
0000006+00  0550756+00 


FDalap  000 

SQRT(V-2+«-2)  0.00086C,  1  000  1  200 

XTtoUaA-  1  Uf*;  ZTioMwk-  1  Unit 
0.000006+00  0572076+00 


FD  Map  1200 

SORT(V~2+H**2)  12.000S6C,  1  600  1  200 

XTfeMufc-  1  Unit;  ZTicMwk-  1  Uni 
0.000006+00  0596766+00 


FDMp  1600 

SQRT(V2+W*2)  16.000S6C,  1  600  1  200 

XTieMaifc-  1  Unk;  Z Tie Maik -  1  Unit 
0.000006+00  0.179296+00 


FDatap  2000 

SORT  (V**2+H**2)  20.000SEC,  1  600  1  200 

X  Tic  Mart*  »  1  Unit;  Z  Tic  Maifc  -  1  Unit 
0.000006+00  0.150706+00 


Rg  (1Hz:2.665km);  Grid  600*500  [dx=0.1km,dt=10ms] 
Profile  15  (seed  91350,  Pr(1)=0.07) 


Figure  16 
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Rg  (1  Hz:2. 665km);  Grid  600*500  [dx=0.1km,dt«10ms] 
Profile  16  (seed  91285.  Pr(1)=0.07) 


SORT  (V*2*+C*2J  Q.00086C,  1  600  1  200 

XTieUwk-  11MI;  ZTteMwk-  1  Ural 


0  000006*00  0  364226*00 


SORT(V2*H~2)  4.000SEC.  1  600  1  200 

X  Tic  Mark  -  1  Unit;  ZTieM»k-  ilk* 


0  000006*00  0.377156*00 


SORT(V**2*H”2)  6.000SEC.  1  600  1  200 

XTfcMaik-  1  U«*;  Z  Tic  Mark-  1  Ur* 


0.000006*00  0264206*00 


FDaUp  1200 

SORT(V~2*fC*2)  12000S6C,  1  600  1  200 

X Tic  Mark-  1  Ur*;  Z  Tie  Mark.  ilk* 
0.000006*00  0.301 356*00 


FD  Map  1600 

SOHX(V~2*H~2)  16.000SEC.  1  600  1  200 

XT*  Mark*  1  Ur*;  Z  Tie  Mark  .  1  Unit 


0.000006*00  0201146*00 


FD  atap  2000 

SORT  (V**2*H**2)  20  000S6C.  1  600  1  200 


X  Tic  Mark-  1  Ur*;  Z  Tic  Mark  -  1  Unit 


0.000006*00  0  194396*00 


Figure  17 
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Rg  (1  Hz:2.665km);  Grid  600*500  [dx-0.1km,dt-10ms] 
Profile  21  (seed  93337,  Pr(1)-0.10) 


SQHT(V*2«M"J)  OOOOS6C.  1  *00  1  MO 
X  Tic  Ma*.  1Ut*.  ZTicMuk.  1  Unit 


0  OOOOOCtOO  0  41700E«00 


FOaMp  400 

SORT(V2««“*)  4  OOOSEC.  1  000  1  000 
XTc  Marti  .  t  lint;  7  Tc  Marti-  1  UM 


0  00000E«O0  0  43740E«00 


SORT(V“2*M**2»  0  OOOSEC,  1  000  1  MO 
XTfeMrak-  1  Ont;  Z  Tic  Marti-  1  UnO 


O.OOOOOEaOO  0.44410E+00 


FD  Map  1200 

80RT<V2«H“2)  120008EC.  1  000  1  200 

X Tic  Maifc*  1  UnO;  ZTic Marti.  1  Ural 


0000006-00  0.361346*00 


FDalap  1000 

SOHT(V-2*H~2)  10.000S6C.  1  000  1  200 


X  Tic  Marti .  1  UnO;  Z  Tic  Marti.  1  Unit 


0  000006*00  0.302326*00 


FDalap  2000 

SORT(V2*H**2)  20  OOOSEC.  1  000  1  200 
X  Tic  Marti  .  1  Un4;  Z  Tie  Maid  -  1  Unit 


0000006*00  0.331646*00 


Figure  18 
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Rg  (1Hz:2.665km);  Grid  600*500  [dx«0.1km,dt»10ms) 
Profile  22  (seed  93347,  Pr(1  )-0.1 0) 


8QRT(v-2-^r*u  aaeMEC,  im  iso 
X  Tic  Marti  -  IMZTcIM-  UM 


0000006*00  0*11  ME >00 


aORT(V-i>M“J)  4.00006C,  1  MO  1  MO 

XTicMnk.  lIMtZTkMnk.  1  Uni 


0.000006*00  040*016*00 


FDaMp  000 

SORTprttH-a  OOOOS6C,  1  OM  1  200 
X Tic Mnk -  1  (Ml;  ZTtoMnk*  1  UnO 


0.000006*00  OM0106*00 


FO  nap  1200 

SORT(V2*H“2)  12.000G6C.  1  000  1  200 

XTIolMi.  1 1MI.  ZTioUnh.  1  Unit 


0000006*00  0.202*6*00 


FDtkp  1000 

90RT(V-J*H-2J  100000EC,  1  000  1  200 

XTkMnk.  UM;  ZTicMafc.  1  Un« 


0.000006*00  0.322076*00 


FOMap  2000 

SQRT(V~2*H~2)  20.000S6C.  1  600  1  200 


XTicMaik.  1  Un»;  Z  Tic  Marti  .  1  Ural 


0.000006*00  0303306*00 


Figure  19 


FD  Map  0 

8QRT(V-a«M-l)  0.0008EC,  1  MO  1  MO 

XTkltek-  1  Uni;  ZTieWwk-  1  Unit 


1 1 1 1 1 1 II I II 1 1 II 1 11 1 1 II 1 1 1 1 1 H 1 1 1 1 II II 1 1 1 1 1 11 1 1 1 1 H 1 1 1 1 1 1 1 1 1 


0  000006+00  0  2100(6+00 


FD  Map  400 

SQRT(V~t*H~2)  4.000S6C.  1  000  1  200 

X  Tc Mwk m  1Unt;  ZTfcUwk-  IUM 
0-000006+00  0  390026+00 


FDMao  000 

8QRT(V"2+H"2)  4.000S6C.  1  too  1  200 

XTkltak-  1 1Mb  ZTIekMo  1  Unit 
0.000006+00  0.3MM6+00 


FDaaap  1000 

8QRT(V“2+W"2)  12.000S6C,  1  (00  1  200 

XTtoUa*.  1  Uni;  ZTicllaak-  1  Uni 
0.000006+00  0  241 206+00 


FD  Map  1000 

SaRT<V2«*r*2)  M.00066C.  1  <00  1  200 

X  Tic  Uvk  •  1  Unit;  ZTb  Mwfc  -  1  Unit 
0.000006+00  02(0556+00 


FD  Map  2000 

SORT(V~2+H"2)  20.000S6C.  1  600  1  200 

X  Tic  Mark-  1  Uni:  Z  Tic  Mark  -  1  Unit 
0.000006+00  0.202426+00 


Rg  (1  Hz:2.665km);  Grid  600*500  {dx«0.1kmtdt«10ms] 
Profile  23  (seed  93356,  Pr(l)-O.IO) 


Figure  20 
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Rg  (1Hz:2.665km);  Grid  600*500  [dx«0.1km,dt-10ms] 
Profile  31  (seed  748546615,  Pr(1)*0.15) 


aORT(V~2+H**2)  OOOStC,  1  MO  1  200 
XTnMUfc-  1  Uni;  Z  Tie  HU*  -  1  Uni 
0.000006*00  0.417896*00 


FDtup  400 

80HT(V~2*W'2)  4.000S6C.  1  000  1  200 

xr«  tkMk.  1  Uni.  Z  Tic  HU*  -  1  Uni 
0000006*00  0  373206*00 


FDUap  MO 

80flr(V2*H~2>  t.OOOSEC,  1  000  1  200 


X  Tie  HU*  -  1  Uni;  Z  Tic  HU*  -  1  Uni 


Ol 000006 *00  020371 E 400 


FDtUp  1200 

9QRT(V**2*H**2)  12000S6C.  1  000  1  200 

XTteMuk-  1  Uni;  Z Tic  HU*-  1  Uni 


0-000006*00  0.101026*00 


FDtUp  1000 

SQRT(V~2*H~2)  16.000S6C.  1  600  1  200 


XTcMwfc-  1  Uni;  Z  Tic  HU*  •  1  Uni 


0.000006*00  0.200566*00 


0M 

wM  ^ 


|  FD  tup  2000 

SORT(V**2*H**2)  20  000SEC.  1  600  1  200 

X  Tic  MaA  -  1  Uni;  Z  Tic  HU*  -  1  Unit 


0.000006*00  0.176206*00 


Figure  21 
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Rg  (1Hz:2.665km);  Grid  600*500  [dx=0.1km,dt-10ms] 
Profile  32  (seed  20333,  Pr(1)«0.15) 


SCWT(V-J*H*-2>  OOOSEC.  t  MO  1 
X  Tic  Mark.  1  Uni;  Z  Tie  Marti  •  1  Unit 


0-000O0E«O0  0.479906 1 00 


FD  Map  400 

SaRT(V2*H“2|  4 OOOSEC,  1  000  1  200 
XTfeMok.  1  Unt;  Z  Tic  Mark.  1 Una 
0000006*00  0.457056*00 


FDtUp  000 

3QHT(V2«H‘*2)  (OOOSEC.  1  000  t  200 


XTieltak.  1  Uni;  Z  Tic  Mark.  1  Ur* 


0000006*00  0200406*00 


ED  (top  1200 

8QRr(V~2*H-2)  12.0006EC.  1  000  1  200 

XTfcMHk-  1  Un«:  Z  Tic  Mark.  1  Unit 


0.000006*00  0202006*00 


ED  (tap  1000 

SQRT(V-2*M~2>  10.000S6C.  1  000  1  200 


XTcMark-  1Un»;  ZTicMaik.  1  Unit 


0.000006*00  0.104116*00 


FO  Map  2000 

SORT<V2*H"2)20.000S6C,  1  600  1  200 


XTic  Mark  -  1  Ur*;  Z  Tic  Mark.  1  Unit 


0.000006*00  0163596*00 


Figure  22 
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Rg  (1Hz:2.665km);  Grid  600*500  [dx*0.1km,dt*10ms] 
Profile  33  (seed  8601 6,  Pr(1  )«0.1 5) 


rUMap  0 

SQRr<V“2*H“2)  0  000SEC.  1  000  1  200 


XTcktek-  IIM-ZTcMmIi.  1  Una 


aoooooe«oo  o.47moe«oo 


FDalap  400 

SORT(V“24*r*2)  4.000S6C.  1  000  1  200 

XTbMmO.  1  Ural.  ZTc  Marti-  1  Ural 
0  000006 *00  0  466656400 


FDalap  000 

SORT(V2*H~2)  6.000SEC,  1  000  1  200 

XTio  Marti  •  1  Ural;  ZTtollwfc-  1  Una 


0.000006 400  0  425706^00 


FDalap  1200 

SORT<V24«~2)  12000SEC,  1  000  1  200 

X  Tie  Marti  -  I  Uni.  Z  Tie  Marti  -  1  Ural 


0.000006400  0295406400 


FDalap  1600 


SORT(V24«"2)  16  000S6C.  1  600  1  200 

XTioMaafc-  1  Ural;  Z Tie  Marti-  1  Unit 


0.000006400  0246226400 


FDalap  2000 

SORT(V~24H**2)  20.000S6C,  1  600  1  200 

XTie  Marti-  1  Un4;  Z Tic  Mark-  1  Unit 


0.000006400  0  279426400 


Figure  23 


0 


FDMap 

SQRT(V**2+W**2)  OOOS6C.  1  100  1  200 
XT* 1  Uni.  Z Tk  Mark  .  1  Out 
0.000006+00  0  36*226  *00 


FD  Map  400 

SORT(V~2+H~2)  4  OOOS6C.  1  600  1  200 
X  Tic  Mark-  lUr*.  ZTkUaA-  1  Out 
0  000006+00  0.330906+00 


FDaMp  000 

SGRT(V"*2+H~2)  O.OOOS6C.  1  000  1  200 
XTtoftlMfc-  1  Unit;  ZTvMmX-  1  Unit 
0.000006+00  0  353076+00 


FDMap  1200 

8QRT(\r*2+H**2)  12  000S6C.  1  600  1  200 

XTioMa*-  1  Unit;  Z Tic Mark -  1  Unit 
0.000006+00  0221406+00 


FDMap  1000 

SQRT(V~2+H^2)  16000S6C,  1  600  1  200 

XTkMa*-  1  Uni;  Z  Tie  Mark-  1  Uni 
0000006+00  0215046+00 


FDMap  2000 

SORT(V*~2+H**2)  20  000S6C.  1  600  1  200 

X  T'c  Mark  -  1  Uni;  Z  Tic  Mark  -  1  Uni 

0.000006+00  0.164366+00 


Rg  (1Hz :2.665km);  Grid  600*500  [dx=0.1km,dt=10ms] 
Profile  42  (seed  93407,  Pr(1)=0.20) 


Figure  25 
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Fomp  o 

SORT<V2**r*2J  o  ooosec.  1  600  1  200 
X Tic  Mark-  1  Unit.  Z  Tic  Mark-  1  Unit 
O.OOOOOE+OO  0.64103E«00 


FOalap  400 

SQRT0T*2*H*‘2)  4.000SEC.  1  600  1  200 

X Tie  Maik-  1  Unit;  Z Tic Maik -  UM 
O.OOOOOE+OO  0.64364E+00 


FOalap  600 

80RT(V**2*H**2)  6.000SEC.  1  600  1  200 

XTkMak-  1  Uni;  Z  Tie  Mark-  1  Unit 
0.00000£«00  0.366S5E«00 


FOalap  1200 

SORT(V**2^H**2)  12.00QSEC,  1  600  1  200 

X  Tie  Maik-  1  Uni;  Z  Tic  Mark-  1  Unit 
O.OOOOOE+OO  0215C2E+00 


FOalap  1600 

SORT(V~24H"2)  16.000SEC.  1  600  1  200 

X  Tic  Maik-  1  Unit;  Z  Tic  Maik  -  1  Unit 
OOOOOOE+OO  0  170996*00 


FOatap  2000 

SORT (V**2+H**2)  20  000SEC.  1  600  1  200 

X  Tic  Mark  -  1  Unit;  Z  Tie  Mark  -  1  Unit 
O  OOOOOE+OO  0.1 6761 E+00 


Rg  (1Hz:2.665km);  Grid  600*500  [dx=0.1km,dt=10ms] 
Profile  43  (seed  12345.  Pr(1)«0.20) 


Figure  26 
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Rg  (1Hz:2.665km);  Grid  600*500  [dx=0.1km,dt«10ms] 
Profile  52  (seed  21696,  Pr(1)=0.25) 


SQRr(V~2«H~2)  ooosec  1  too  1  *00 
XTcUa*«  1  On*.  Z Tie kbtfk  .  I  UM 
0  000006*00  0.47SNE«00 


FD  Map  400 

SOfn(V**2»H"2)  4  OOOSEC  1  000  1  200 
XTte  Marie  1Un»;  Z  Tic  Marie  1  Una 
0  000006*00  0  543286  *00 


FOatap  000 

SORT(V*2*H“2)  8  000S6C  1  «00  1  200 


XTkMarte  lUnt:  ZTkMaik.  t  Unit 


0  000006*00  0242446*00 


FDtUp  1200 

8QRT(V“2*H"2)  120008EC.  1  MO  1  200 

X  Tic  Marie  1  Uni;  Z  Tic  Marie  1  Unit 
0  000006*00  01 75446*00 


FDatap  1800 

SOFTT(V~2*Fr*Z)  18.000SEC.  1  too  1  200 

X  Tie  Marie  1  Uni;  ZTic  Mark  .  1  Unit 


O.OOOOOETOO  0.20668E+00 


FD  Stop  2000 

SQRT(V**2tH**2)  20. OOOSEC,  1  000  1  200 

XTicMirtc  1  Unt;  ZTic  Mark-  1  Unit 


O.OOOOOEtOO  067041 E -01 


Figure  28 


FOfltop  0 

SQRr(V-l*H~J)  00066C.  I  «00  I  MO 
XTcUak-  1  Unt;  Z  Tie  Mast  •  1  UM 
0000006400  0.551236400 


FO«Mp  400 

90RT(V-24W2)  4.00086C.  1  800  1  200 

XTcMwk-  1  Unt;  ZTteMaik-  1  Uni 
0  000006 400  0.580506400 


FDMap  MO 

acwr(v-24H"2)  IOOOS6C.  i  too  i  zoo 
XTtoMnk-  lUnt;  ZTfelMi-  1  Uni 
0000006 400  0347016400 


FOaUp  1200 

90RT(V~24«"2)  12000S6C.  1  800  1  ZOO 

X Tio t4nk -  lUnt;  ZTlcMnk-  1  Unit 
0.300006400  0301846400 


FO«Up  1000 

SQHT(V~Z4»r*Z)  18000S6C.  1  800  1  zoo 

XTeMnk-  1  Uflt;  ZTicMnV-  1  Unit 
0  000006 400  0.217796400 


FD  slap  2000 

SQRT(V-Z4H"Z)  20  000SEC.  1  800  1  200 

X  Tic  Mark  *  1  Ur*.  Z  Te  Mark  .  1  Unit 
0.000006400  0  992036-01 


Rg  (1Hz:2.665km);  Grid  600*500  [dx=0.1km,dt.=1  Oms] 
Profile  53  (seed  86037,  Pr(1)=0.25) 


Figure  29 
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Rg  (1  Hz:2.665km);  Grid  600*500  [dx«0.1km,dt-10ms] 
Profile  61  (seed  93459,  Pr(1)-0.30) 


GQRT(V“2*H~25  00066C.  1  *00  1 

XTtoMak.  1  Uni;  Z Tie MaA •  1  Uni 


o.oooooe«oo  o  MioaEtoo 


SQRT(V~2*H~2)  4  000 SEC,  1  >00  1  200 

X  re  IM  ■  1  Unit;  2  Tic  Mart-  1  Unit 


0000006*00  0655526*00 


FOaMp  tOO 

80RT(V~2*H"2)  0.00086C,  1  600  1  200 

XTkMMk-  1  Uni.  ZTIcMaik.  Unit 


0000006*00  0.000066*00 


FD  Map  1200 

aafrr(v~2+H~2)  naxeec.  i  eoo  i  200 

XTIoMMk-  1  Ur*:  ZTtolvUik.  1  Unit 


0.000006*00  0.400726*00 


FDalap  1000 

SORT<V2*H~2)  10.000S6C.  1  000  1  200 

XTfcMwk-  1  Un4:  ZTIcMaik.  1  Unit 


0.000006*00  0.420636*00 


FDatap  2000 

SQRT(V~2*ft~2)  20.000S6C,  1  600  1  200 

XTicMaik-  1  Unit:  Z  Tic  Marie-  1  Unit 
0.000006*00  0.106356*00 


Figure  30 
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FDatap  0 

SOHT<V~2+H~2»  OOOSEC.  1  900  1  200 
X  Tic  Marti  -  1  Un*.  Z  Tic  Merit-  1  Ural 
0  000006*00  0  941036+00 


FO  Map  400 

S0RT<V2+H**2>  4  OOOSEC.  1  800  1  200 
X Tic  Marti  -  1  Unit;  ZTie  Merit*  1  Ural 
0.000006+00  0  802126+00 


FD  Map  800 

SQRT(V~2+H*a2)  9  OOOSEC.  1  800  1  200 
X  Tie  Merit*  1  Ural;  ZTieMerit-  1  Unit 
0.000006+00  0  448586+00 


FDmp  1200 

SQRT<V2+H~2)  12  OOOSEC.  1  800  1  200 
XTeMek.  1  Unit;  Z  Tie  Merit-  1  Ural 
0000006+00  0.3061 7E +00 


FD  step  1800 

SQRT(V**2+H**2)  18. OOOSEC.  1  800  1  200 

X  Tie  Merit*  1  Unit;  Z  Tic  Merit  -  1  Unit 
0  000006+00  0  317796+00 


FD  step  2000 

SOnT(V~2+H~2)  20  000S6C.  1  600  1  200 

X  Tic  Merit  *  1  Unit;  Z  Tic  Merit  *  1  Unit 
0  000006+00  0  123306+00 


Rg  (1Hz:2.665km);  Grid  600*500  [dx=0.1km,d1=10ms] 
Profile  62  (seed  93468,  Pr(1)=0.30) 


Figure  31 


•a«T(v-7»*r-2)  aoooesc.  i  aoe  i 

XTaMa*.  1  Uni;  Z  Tto  DM  -  t  UM 
aoooooe*oo  0  041836.00 


aaRT(v-j»tn)  4  000 see.  i  too  i 

XTkUak-  1 UM.  Z  Tic  MM-  iUm 
0000006*01  0I7WH  lOO 


SORT (V-l^TT)  IOOOOCC.  1  MO  1  100 
XTtoMMi-  UMtZTkM**-  1  UM 
0.000006*410  O  SM  176*00 
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FOMp  1100 

8QRT(V~l4H~a)  1200006C.  1  MO  1  200 

X Tic IM -  IUM;  ZTbMa*.  IUM 
0000006*00  02SS06E«00 
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FOOo  1000 

8QHT(V-2^r-2)  1000066C,  1  000  1  200 

XTfcMvk-  1  UM;  ZTloMuk-  IUM 
0000006*00  0190576*00 
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iFDMp  2000 

SOflT(V~2*H~2) 200008CC.  1  600  1  200 

XTicMiA-  1  UM.  ZTieMe*  -  IUM 


Rg  {1Hz2.665km);  Grid  600*500  [dx«0.1km,dt=10ms] 
Profile  63  (seed  93479,  Pr(1)«0.30) 


Figure  32 
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Rg  (1Hz2.665km);  Grid  600*500  [dx-0.1km.dt-  10ms] 
Profile  64  (seed  55473,  Pr(1)-0.32) 


SQRT(V-2*H“Q  OOOS6C.  1  *00  I  MO 
XTieMa*-  1  OM.  Z  Te  Marti  .  IUM 


0000006-00  0  041036.00 


FDaNp  400 

SORT <V2-H“2)  4.  OOOSEC,  1  000  1  200 

X  Tic  Marti-  1  UnO.  ZTc  Marti  -  1  Unit 


0  000004*00  0  724006-00 


FOalap  000 

SQHT(V-*«M“2)  0.000S6C,  1  000  1  200 

XTleMMk-  1  UnO;  Z  Tic  Marti  -  1  Unit 


0000006-00  0.310206-00 


FD  nap  1200 

80Rr(V-2*H~2)  12.000S6C.  1  000  1  200 

X  Tic  Marti-  1 1M.  zr- Marti-  1  Una 


O.OOOOOE«00  0 20201 E -00 


FDatap  1000 

9QRT(V~2-W”2)  1000086C.  1  000  1  200 

X  Tie  Marti-  1  Un«:  Z  Tc  Marti .  1  Una 


0  000006-00  0205576.00 


FDatap  2000 

SQRT(V~2-«-*21 20000SEC,  1  000  1  200 


X Tic  Marti*  1  Unt;  ZTic Marti  •  1  Una 


0  000006-00  0  51740601 


Figure  33 
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SQRT(\T*2+H**2)  OOOSEC.  1  000  1  200 

X Telia*-  1  Uni;  Z Tc Ma* -  1  Ural 
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FOatap  0 

8QRT(V**«H"*}  OOOSEC.  1  *00  <  200 
XTtoltak.  1  UM.  2 Tic  Uwk  .  1IM 
0.000006400  0.631566*00 


FD  Map  400 

8QRT(V“2*H-2)  4.000SEC,  1  600  1  200 

XTieMMk-  lOni;  ZTcMaik.  1UM 
0000006*00  0.536476*00 


FD  aup  600 

aORT(V“2*+r-2)  6  OOOSEC.  1  600  1  200 
XTieMMk-  1  Link;  ZTkMwk.  1  Unit 
0.000006*00  0301666*00 


FDWp  1200 

8Qf«rpn*H~2)  12  OOOSEC.  1  600  1  200 
XTbtlaik.  1  Uni.  Z Tc Maik .  1  Unit 
0000006*00  0222*1 E *00 


FDalap  1600 

SORT(V~2*H"2)  16000SEC,  1  600  1  200 

X  Tic  Mai*.  1  Uni:  ZTteMaik-  1  Unit 
0.000006*00  0.212106*00 


FOatap  2000 

SORT(V2*H"2)  20  OOOSEC,  1  600  1  200 
XTcMartc  -  1  On*:  ZTcMaifc*  1  Ural 
0.000006*00  0.666146-01 


Rg  (1  Hz:2.665km);  Grid  600*500  [dx-O.lkm.dt-IOms] 
Profile  71  (seed  22192,  Pr(1)=0.35) 


Figure  36 
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FD  Map  0 

SQRt(v-*^r-a)  ooosec.  i  ooo  i  zoo 

XTieklMk-  1  Una.  Z Tn Mat  -  IUM 
OOOOOOEtOO  0  417006*00 


FOaMp  400 

SQRT(Vl4H“2>  » OOOS6C.  1  000  1  200 
XTteMnk.  1  Ural.  ZTcMat-  1  Uni 
O.OOOOOEtOO  0  370306*00 


FOaMp  *00 

8QRT(V“2*H"2)  0.  OOOSEC,  1  000  1  200 

XTieMnk.  1  UM;  Z Tie Mat .  1  Ural 
0  00000E*O0  01 50576*00 


FOaMp  1200 

SORT(V-2*H-2)  12000S6C.  1  000  1  200 

XTieMnk-  1  Un*.  ZTtaMaik.  1  Unt 
0000006*00  0.110426*00 


FOaMp  1000 

SQRT(\r*2*H*'2)  10.0006EC.  1  000  1  200 

XTieMnk-  1  Uni;  ZTioMnk.  1  Ural 
O.OOOOQE«00  0.117016*00 


FDatap  2000 

SORT  (V*2*H"2)  20. OOOSEC.  1  000  1  200 

XTieMnk-  1  Unk;  ZTicMuk-  1  Unit 
0.000006+00  0  48721E-01 


Rg  (1Hz:2.665km);  Grid  600*500  (dx=0.1km,dt-10ms] 
Profile  72  (seed  31087,  Pr(1)=0.35) 


Figure  37 


FD*ap  0 


min . . 


SQffr<V~2*H~2)  0  OOOSEC.  1  COO  1  200 
XTb  Muk  •  1  Ural,  Z  Tic  Umtk  -  1  Ural 
0.000006400  0-304226*00 


FDMp  400 

SORT(V"*2*W~2)  4  OOOSEC.  1  600  1  200 
XTcMwk-  1  Uni;  Z  Tie  Mafk  •  1  Ural 
0000006*00  0.510256*00 


FO  Mtp  600 

SQRT(V*2+HM2)  6  OOOSEC.  1  600  1  200 
XTfcMrak-  1  Uni;  ZTtoMufc*  1  Unit 
0000006*00  0.166666*00 


FO  atop  1200 

SORT(V2**r*2>  12  OOOSEC.  1  600  1  200 
X  Tic  Marti  ■  1  Untt:  ZTioMrak-  1  Unit 
0.000006*00  0 105576*00 


FDstap  1600 

SORT(V~2*H"2)  16  OOOSEC,  1  600  1  200 
XTfcMwk.  1  Ur*;  Z  Tie  U*k  -  1  Unit 
0.000006*00  0.51 5996-01 


FOstap  2000 

SORT(V2*H**2)  20  OOOSEC.  1  600  1  200 
XTic  Maifc  •  1  Ural;  Z  Tic  M*ik-  1  Unit 
0  000006*00  0  338746-01 


Rg  (1Hz:2.665km);  Grid  600*500  [dx-0.1km.dta  10ms] 
Profile  73  (seed  86061  ,  Pr(1)-0.35) 


Figure  38 
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le  11  (seed  22887,  Pr(1)=0.05),  Surface  Sensors 

same  scale,  decimated  by  S 


0.162E+00 

0.156E+00 


le  13(seed  85995,  Pr(1)=0.05),  Surface  Sensors 

same  scale,  decimated  by  5 


le  14(seed  91333,  Pr(1)=0.07),  Surface  Sensors 
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Profile  15(seed  91350,  Pr(1)=0.07),  Surface  Sensors 

VERT,  same  scale,  decimated  by  5 
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Profile  21  (seed  93337,  Pr(1)=0.10),  Surface  Sensors 
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0.197E+00 


sd  93347,  Pr(1)=0.10),  Surface  Sensors 
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8d  20333,  Pr(1)=0.15),  Surface  Sensors 
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Profile  33(seed  86016,  Pr(1)=0.15),  Surface  Sensors 


0.181E+00 

0.163E+00 

0.163E+00 


Profile  41  (seed  93397,  Pr(1)=0.20),  Surface  Sensors 
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00+30910 
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Profile  43(seed  12345,  Pr(1)=0.20),  Surface  Sensors 

VERT,  samo  scale,  docimaled  by  5 


0.166E+00 

0.173E+00 


sd  21696,  Pr(1)=0.25),  Surface  Sensors 


0.156E+00 

0.160E+00 


ad  86037,  Pr(1)=0.25),  Surface  Sensors 


0.81  #E -01 


Bd  93459,  Pr(1)=0.30),  Surface  Sensors 


00+3«SVO 
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sd  93468,  Pr(1)=0.30),  Surface  Sensors 


0.112E+00 

0.112E+00 


Profile  63(seed  93479,  Pr(1)=0.30),  Surface  Sensors 


0.105E+00 

0.855E-01 


VERT,  same  scale,  decimated  by  5 


0.161E+00 

0.144E+00 


0.174E+00 

0.157E+00 


0O+3W10 

00*31610 


Profile  72(seed  31087,  Pr(1)=0.35),  Surface  Sensors 


Peak  Z-0.88670  ;  Peak  H-0.71420  ;  Scaled 


Particle  Motion,  Trace  24  Model  00,  Rg  (1  Hz:  2.665km/s) 


Particle  Motion,  Trace  59  Model  33,  Scattered  Synthetics  (Rotated) 


0.0  0.0  5.0  10.0  15.0  20.0 

Particle  Motion,  Trace  59  Model  72,  Scattered  Synthetics  (Rotated) 


Variation  vs.  Gamma  Ll-Varlatlon  vs.  Transmitted  Power 


Standard  LS:  S*-0.393(0.056),  1*0.280(0.020),  rho*-0.81  Standard  LS:  8*0.894(0.1 04),  1*0.41 3(0.036),  rtw*047 


L2-Variation  vs.  Gamma  L2-Varlatlon  vs.  Transmitted  Power 
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Standard  US:  S«-7.205(1.118),  l<04S4(0.037),  rho«-0.7»  Standard  LS:  S-17.118(1.e71),  1.0.154(0.061),  rtnnO.ee 


RMS  Elevation  vs.  Gamma  RMS  Elevation  vs.  Transmitted  Power 


ZH1  teeuiuieg 


juapujaoo  uoipeuau 


Standard  LS:  S--0.1 27(0.045),  1.0.233(0.030),  rho»-0.49  Standard  LS:  8.0.326(0.089),  U0 .490(0.080),  rt*o.O M 


Durbln-Watson  Number  vs.  Gamma  Durbln-Watson  Number  vs.  Transmitted  Power 


Standard  LS:  S»1  -876(1.527),  1.0.132(0.024),  rho.0.24  Standard  LS:  8.-4.545(3.236).  1-0.754(0.052),  rho-0.27 


Peak  Wavenumber  vs.  Gamma  Peak  Wavenumber  vs.  Transmitted  Power 
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Standard  LS:  S--0.1 38(0.032),  1.0.248(0.024),  rtio«-0.65  Standard  LS:  S^)  .263(0.072),  U0.51 8(0.055),  rtto-0  M 


Corners  per  KM  vs.  Gamma  Comers  per  KM  vs.  Transmitted  Power 


Corners  per  KM  Comers  per  KM 

Standard  LS:  S--0.070(0.020),  1*0.301(0.044),  rtio-0.57  Standard  LS:  8*0.202(0.035),  1*0.281(0.074),  rho-0.76 


Normalized  Arc  Length  vs.  Gamma  Normalized  Arc  Length  vs.  Tra 


Standard  LS:  S— 0.949(0.136),  1.1.229(0.154),  rtw>.-0.81  Standard  LS:  S-2.1 59(0.251 >,  1-1.746(0.2*4),  ihOmOM 


A  Priori  Pr[Elevatlon  Change]  vs.  Gamma  A  Priori  Pr[Elevation  Change]  vs.  T 
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Standard  LS:  S»-0.566<0.085), 1-0.269(0.019).  rho-0.80  Standard  LS:  8-1^00(0.156),  1-0.437(0.015),  rtto-0.66 


P,  Model  11  (seed  22887,  Pr(l)=o.05) 

vert,  same  scale,  decimated  by  5 
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S.  Model  11  (seed  22887,  Pr(1)=0.05) 

hori,  same  scale,  decimated  by  5 


40.0  Second 


Figure  103 
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Second 


P,  Model  12(seed  22C92,  Pr(1)-0.05) 

vert,  same  scale,  decimated  by  S 
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Second 


S,  Model  12(seed  22892,  Pr(1)=0.05) 

hori,  same  scale,  decimated  by  5 


Figure  104 
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P,  Model  13(seed  85995,  Pr(1)=0.05) 

vert,  same  scale,  decimated  by  5 
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0.0  10.0  20.0 
S.  Model  13(seed  85995,  Pr(1)=0.05) 

hort,  same  scale,  decimated  by  5 
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Figure  105 
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P,  Model  14(seed  91333,  Pr(1)*0.07) 


vert,  same  scale,  decimated  by  5 
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S.  Model  14  (seed  91333,  Pr(i)=0.07) 


hori,  same  scale,  decimated  by  5 


Figure  106 
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P.  Model  16(seed  91285,  Pr(1)«0.07) 

vert,  same  scale,  decimated  by  5 
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S,  Model  16(seed  91285,  Pr(1)=0.07) 


hori,  same  scale,  decimated  by  5 


Figure  108 
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P,  Model  21  (seed  93337,  Pr(1)»0.10) 

vert,  same  scale,  decimated  by  5 


0.0  10.0 
S,  Model  21  (seed  93337,  Pr(1)=0.10) 

hori,  same  scale,  decimated  by  5 
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Figure  109 
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0.0  10.0  20.0  30.0  40.0 

P.  Model  22(seed  93347.  Pr(l)-O.IO) 

vert,  same  scale,  decimated  by  5 


S.  Model  22{seed  93347.  Pr(1)=0.10) 


hori,  same  scale,  decimated  by  5 
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Second 


Figure  110 
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P,  Model  23 (seed  93356.  Pr(l)-O.IO) 

vert,  same  scale,  decimated  by  5 
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Second 


hori,  same  scale,  decimated  by  5 


Figure  111 
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31 


P.  Model  31  (seed  748546615.  Pr(1)«0.15) 
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14.1  0.148E-01 

17.1  0.160E-01 

19.9  0.172E-01 

22.6  0.184E-01 

25.3  0196E-01 

27.8  0.207E-01 

30.3  0.217E-01 

32.7  0.227E-01 

35.0  0.238E-01 

36.9  0.247E-01 

38.3  0.256E-01 

40.1  0.264E-01 

41.5  0S6BE-01 

42.9  0S72E-01 

44S  0S74E-01 

45.1  0S78E-01 

54.8  0.228E-01 


Second 


vert,  same  scale,  decimated  by  5 


S,  Model  31  (seed  748546615.  Pr(1)=0.15) 


-19.6  0.414E-01 

-20.7  0.409E-01 

-21.3  0.444E-01 

-15.7  0.412E-01 

-13.8  0.434E-01 

-9.7  0.461 E-01 

-7.7  0.465E-01 

-5.1  0.490E-01 

-1.4  0.542E-01 

-OS  0.534E-01 

5.0  0.543E-01 

8.0  0.543E-01 

11.2  0.542E-01 

14.1  0.514E-01 

17.1  0.493E-01 

19.9  0.475E-01 

22.6  0.461 E-01 

25.3  0.396E-01 

27.8  0.406E-01 

30.3  0.406E-01 

32.7  0.308E-01 

35.0  0.305E-01 

36.9  0.259E-01 

38.3  0.342E-01 

40.1  0.463E-01 

41.5  0.473E-01 

42.9  0.371  E-01 

44.2  0.241  E-01 

45.1  0.1 61  E-01 

54.8  0.221  E-01 


Second 


hori,  same  scale,  decimated  by  5 


Figure  112 
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P.  Model  32(seed  20333,  Pr(1)»0.l5) 

vert,  same  scale,  decimated  by  5 


0.226E -01 
0.230E-01 
0.233E-01 
0.225E-01 
0  223b -01 
0.21 4E -01 
0.209E-01 
0.201  E-01 
0.191E-01 
0  181  E-01 
0.170E-01 
0.150E-01 
0.148E-01 
0.142E-01 
0.126E-01 
0.116E-01 
0.133E-01 
0.122E-01 
0.158E-01 
0.140E-01 
0.132E-01 
0.125E-01 
0.140E-01 
0.123E-01 
0.138E-01 
0.120E-01 
0.128E-01 
0.168E-01 
0.132E-01 
0.100E-01 


40.0  Second 


0.334E-01 
0.345E-01 
0.382E-01 
0.355E-01 
0.422E-01 
0.504E-01 
0.577E-01 
0.559E-01 
0.570E-01 
0.564E-01 
0.571  E-01 
0.551  E-01 
0.540E-01 
0.551  E-01 
0.529E-01 
0.501  E-01 
0.499E-01 
0.456E-01 
0.348E-01 
0.361  E-01 
0.354E-01 
0.356E-01 
0.416E-01 
0.448E-01 
0.456E-01 
0.421  E-01 
0.390E-01 
0.331  E-01 
0.249E-01 
0.207E-01 


S,  Model  32(seed  20333,  Pr(1)=0.15) 

hori,  same  scale,  decimated  by  5 


40.0  Second 


Figure  113 
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vert,  same  scale,  decimated  by  5 

J-L  1  I  -L  I  I  I  I  1  I  I  I  I  I  1  1  I  I  1  I  I 


I  II  I  I  I  I  II  |  I  I  I  I  I  I  I  I  I  |  I  | 
>0  10.0  20.0 
S.  Model  33 (seed  86016.  Pr(1)-0.15) 

hori,  same  scale,  decimated  by  5 


•23.6 

0.1676-01 

-196 

0.151E-01 

-20.4 

0.1216-01 

-156 

0.120E-01 

-115 

0.117E-01 

-3.5 

0.108E-01 

-5.9 

0.9906-02 

-22 

0.939E-02 

1.9 

0.9326-02 

5.8 

0.132E-01 

9.6 

0.1256-01 

125 

0.1476-01 

14.6 

0.1356-01 

17.7 

0.1486-01 

202 

0.16SE-01 

22.8 

0.1506-01 

252 

0.1346-01 

275 

0.1926-01 

29.6 

0.1346-01 

315 

0.1556-01 

33.4 

02006-01 

355 

0.1576-01 

37.6 

0.162E-01 

395 

0.1666-01 

41.5 

0.1696-01 

432 

0.171E-01 

44.8 

0.1716-01 

45.1 

0.171E-01 

452 

0.1806-01 

55.4 

0.1406-01 

Second 


-23.6 

0.379E-01 

-195 

0.3706-01 

-20.4 

0.380E-01 

-156 

0.3916-01 

-115 

0.3616-01 

-8.5 

0.3776-01 

-5.9 

0.4496-01 

-22 

0.4906-01 

1.9 

0.5206-01 

56 

05726-01 

9.6 

0.5846-01 

12.5 

0.6126-01 

146 

0.6476-01 

17.7 

0.6336-01 

202 

0.578E-01 

22.8 

0.5286-01 

252 

0.4976-01 

275 

0.4976-01 

296 

0.5016-01 

31.5 

0.5166-01 

33.4 

0.5076-01 

355 

0.4836-01 

37.6 

0.4576-01 

395 

0.3806-01 

41.5 

0.327E-01 

432 

0.310E-01 

44.8 

0.262E-01 

45.1 

0.345E-01 

45.2 

0.402E-01 

55.4 

0.429E-01 

30.0 


40.0 


Second 


Figure  114 


0.1146-01 
0.9886-02 
0.1496-01 
0.102E-01 
0.8656-02 
0.7886-02 
0  7936-02 
0.8996-02 
0.1096-01 
0.9686-02 
0.1086-01 
0.1286-01 
0.1196-01 
0.1756-01 
02206-01 
0.1796-01 
0.1926-01 
0.1466-01 
0.1536-01 
02016-01 
02046-01 
0.1746-01 
0.1706-01 
02146-01 
0.1826-01 
0.1956-01 
0.1926-01 
02196-01 
02336-01 
0.1756-01 


P,  Model  41  (seed  93397,  Pr(l)-0.20) 

vart,  aarne  seal*,  decimatad  by  5 


40.0  Second 


0.5356-01 

0.4626-01 

0.4566-01 

0.4616-01 

0.4556-01 

0.4746-01 

0.4966-01 

0.5016-01 

0.5286-01 

0.5216-01 

0.5276-01 

0.5236-01 

0.4966-01 

0.4746-01 

0.562E-01 

0.6716-01 

0.5446-01 

0.410E-01 

0.3486-01 

0.3686-01 

0.3306-01 

0.2996-01 

0.3196-01 

0.3106-01 

0.3306-01 

0.3476-01 

0.4146-01 

0.320E-01 

0.3576-01 

0.342E-01 


S,  Model  41  (seed  93397,  Pr(1)=0.20) 

hori,  same  scale,  decimated  by  5 


40.0  Second 


Figure  115 
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vert,  same  scale,  decimated  by  5 


I  I  I  I  I  I  I  I  |  I  I  I  I  I  I  I 
0.0  10.0 
S,  Model  42(seed  93407,  Pr(1)=0.20) 

hori,  sane  scale,  decimated  by  5 


-23.6 
-17.7 
-18.3 
-17.7 
-13.1 
-7.8 
-6.2 
-2.6 
0.4 
3.5 
7  5 
8.0 
11.9 

16.4 
18.8 
21.2 

23.4 

25.6 

27.8 
30.0 

32.1 

34.1 
36.0 

37.9 
39.8 

41.6 

43.4 
45.0 

47.2 

55.2 


0.171E-01 
0.138E-01 
0.158E-01 
0.181E-01 
0  1 696-01 
0.161E-01 
0  1 626 -01 
0.159E-01 
0.135E-01 
0.138E-01 
0.138E-01 
0.132E-01 
0.149E-01 
0.168E-01 
0.185E-01 
0501 E-01 
0  216E-01 
0.229E-01 
0  241  E-01 
0551  E-01 
0.257E-01 
0 .2606-01 
0  2626 -01 
0561E-01 
05S9E-01 
0 .2636-01 
0.2666-01 
05676-01 
0565E-01 
0518E-01 


Second 


-23.6 

0.386E-01 

-17.7 

0.395E-01 

-185 

0.417E-01 

-17.7 

0.493E-01 

-13.1 

0.521  E-01 

-7.8 

0.550E-01 

-65 

0.575E-01 

-2.6 

0.567E-01 

0.4 

0.6156-01 

3.5 

0.6246-01 

75 

0.680E-01 

8.9 

0.691  E-01 

11.9 

0.685E-01 

16.4 

0.654E-01 

18.8 

0.657E-01 

21.2 

0.5386-01 

23.4 

0.518E-01 

25.6 

0.583E-01 

27.8 

0.500E-01 

30.0 

0.442E-01 

32.1 

0.418E-01 

34.1 

0.347E-01 

36.0 

0.317E-01 

37.9 

0.341  E-01 

39.8 

0.350E-01 

41.6 

0.502E-01 

43.4 

0.556E-01 

45.0 

0.514E-01 

47.2 

0.431  E-01 

55.2 

0.362E-01 

20.0 


30.0 


40.0 


Second 


Figure  116 


1.0  10.0 
P.  Model  51  (seed  21691,  Pr(l)«0.25) 

vert,  same  scale,  decimated  by  5 


0.0  10.0 
S,  Model  51  (seed  21691 ,  Pr(l)=0.25) 

hori,  same  scale,  decimated  by  5 


0 ,232b -Ot 
0.239E-01 
0.237E-01 
0.239E-01 
0.233E-01 
0.227E-01 
0.221  E-Ot 
0.21  IE-01 
0.200E-01 
0.189E-01 
0.180E-01 
0.170E-01 
0.163E-01 
0.159E-01 
0.161E-01 
0.163E-01 
0.166E-01 
0.168E-01 
0.171E-01 
0.173E-01 
0  174E-01 
0.175E-01 
0.175E-01 
0.175E-01 
0.174E-01 
0.174E-01 
0.173E-01 
0.173E-01 
0.174E-01 
0.145E-01 


40.0  Second 


0.358E-01 
0.383E-01 
0.405E-01 
0.453E-01 
0.517E-01 
0.570E-01 
0.61  IE-01 
0.608E-01 
0.57  IE-01 
0.571E-01 
0.533E-01 
0.463E-01 
0.371  E-01 
0.327E-01 
0.366E-01 
0.284E-01 
0.257E-01 
0.307E-01 
0  290E-01 
0.299E-01 
0.360E-01 
0.354E-01 
0.310E-01 
0.329E-01 
0.385E-01 
0.382E-01 
0.398E-01 
0.292E-01 
0256E-01 
0.246E-01 


40.0  Second 


Figure  118 


P,  Model  52(seed  21696,  Pr(l)=0.25) 

vert,  same  scale,  decimated  by  5 


0  182E-01 
0.159E-01 
0.163E-01 
0.161E-01 
0.173E-01 
0.158E-01 
0154E-01 
0.144E-01 
0.158E-01 
0.173E-01 
0  183E-01 
0190E-01 
0194E-01 
0.194E-01 
0.192E-01 
0.189E-0I 
0.185E-01 
0.182E-01 
0.180E  Ot 
0.176E  01 
0.171E-01 
0  169E-01 
0.168E-01 
0.173E-01 
0.180E-01 
0.189E-01 
0.198E-01 
0.204E-01 
0.221  E-01 
0.175E-01 


40.0  Second 


0.378E-01 

0.332E-01 

0.372E-01 

0.316E-01 

0.349E-01 

0.371  E-01 

0.401E-01 

0.407E-01 

0.364E01 

0.334E-01 

0.465E-01 

0.394E-01 

0.429E-01 

0.501  E-01 

0.517E-01 

0.547E-01 

0.593E-01 

0.626E-01 

0.648E-01 

0.652E-01 

0.646E-01 

0.580E-01 

0.494E-01 

0.433E-01 

0.406E-01 

0.375E-01 

0.340E-01 

0.287E-01 

0.287E-01 

0.213E-01 


S,  Model  52(seed  21696,  Pr(1)=0.25) 

hori,  same  scale,  decimated  by  5 


Figure  119 
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-23.1 

0.243E-01 

-22.6 

0.189E-01 

-20.9 

0.174E-01 

-16.3 

0.171E-01 

-15.1 

0.169E-01 

-11.1 

0  179E-01 

-7.1 

0.179E-01 

-5.0 

0.176E-01 

-1.5 

0.171E-01 

2.7 

0.166E-01 

6.1 

0.160E-01 

9.6 

0.152E-01 

13.0 

0.153E-01 

17.1 

0.164E-01 

20.4 

0.171E-01 

23.3 

0.143E-01 

25.6 

0.145E-01 

27.5 

0.140E-01 

29.3 

0.150E-01 

31.2 

0.138E-01 

33.1 

0.149E-01 

34.4 

0.166E-01 

35.3 

0.186E-01 

37.4 

0.204E-01 

38.7 

0.220E-01 

40.4 

0.230E-01 

41.8 

0.241  E-01 

44.6 

0.236E-01 

43.9 

0.258E-01 

54.7 

0.198E-01 

Second 


P,  Model  53(seed  86037,  Pr(1)=0.25) 


vert  same  scale,  decimated  by  S 


-23.1 

0.661  E-01 

-22.6 

0.589E-01 

-20.9 

0.591  E-01 

-16.3 

0.582E-01 

-15.1 

0.561  E-01 

-11.1 

0.547E-01 

-7.1 

0.486E-01 

-5.0 

0.450E-01 

-1.5 

0.431E-01 

2.7 

0.444E-01 

6.1 

0.416E-01 

9.6 

0.427E-01 

13.0 

0.438E-01 

17.1 

0.358E-01 

20.4 

0.331  E-01 

23.3 

0.359E-01 

25.6 

0.361  E-01 

27.5 

0.365E-01 

29.3 

0.443E-01 

31.2 

0.404E-01 

33.1 

0.357E-01 

34.4 

0.414E-01 

35.3 

0.407E-01 

37.4 

0.480E-01 

38.7 

0.517E-01 

40.4 

0.503E-01 

41.8 

0.452E-01 

44.6 

0.302E-01 

43.9 

0.237E-01 

54.7 

0.286E-01 

Second 


S,  Model  53(seed  86037,  Pr(1)=0.25) 


hori,  same  scale,  decimated  by  5 


r 


Figure  120 
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-22.6 

0J78E-O1 

-18.8 

0.220E-01 

-20.0 

0.271E-01 

-145 

0.239E-01 

-9.8 

0212E-01 

-7  0 

0.212E-01 

5.7 

0.214E-01 

-7.4 
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-0.2 

0.195E-01 

3.4 

0.148E-01 
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0.145E-01 
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0.160E-01 

13.7 

0.169E-01 

16.8 
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19.8 

0.1B1E-01 
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0.194E-01 

25.3 
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0.257E-01 
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0.217E-01 

33.1 

0.265E-01 

35.4 
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0.225E-01 

39.8 

0.217E-01 

41.9 

0.257E-01 

44.0 

0.363E-01 

45.8 

0.360E-01 

47.4 

0275E-01 

48.5 

0.270E-01 

49.9 

0.234E-01 

58.1 

0.170E-01 

0.0 


10.0 


20.0 


30.0 


40.0  Second 


P,  Model  61  (seed  93459.  Pr(1)-0.30) 


31 

32 

33 

34 
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52 
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57 

58 

59 

60 


vert,  same  scale,  decimated  by  5 

I  I  I  1  1  I  I  I  I  1  I  1  I  I  M  I  l  I  1 


0.0 


I  I  I  I  I  I  I  I  I  |  I  I  I  I  I  NTT 

10.0  20.0 


S.  Model  61  (seed  93459,  Pr(1)=0.30) 

hori,  same  scale,  decimated  by  5 


-22.6 

0.380E-01 

-18.8 

0.445E-01 

-20.0 

0.512E-01 

-145 

0.581  E-01 

-9.8 

0.593E-01 

-7.0 

0.622E-01 

-5.7 
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-7.4 

0.546E-01 
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0.489E-01 
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0.466E-01 

41.9 
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0.486E-01 

i  i  i  i  |  i  i  n  ii  i  i  i | 

30.0  40.0 


Second 


Figure  121 
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Second 


P.  Model  62(seed  93468,  Pr(1)«0.30) 

vert,  same  scale,  decimated  by  5 
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P.  Model  72(seed  31087,  Pr(l)-0.35) 

vert,  same  scale,  decimated  by  5 
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P.  Model  73(seed  86061 .  Pr(l)«0.35) 
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